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— ^  A  study  of  the  wave  induced  springing  of  the  Great  Lakes  ore  carrier  M/V  STEWART 
J.  CORT  is  being  done  in  two  phases.  Phase  one  is  the  installation,  test,  and  evalua¬ 
tion  of  a  shipboard  wave  measurement  system.  Phase  two  is  the  measurement  and 
correlation  of  the  waves  encountered  by  the  vessel  and  the  vessel's  response.  This 
report  describes  the  wave  measurement  system  and  presents  some  of  the  test  results. 
The  wave  measurement  system  consists  of  two  radars  used  to  profile  the  surface  waves 
and  two  sets  of  accelerometers  to  record  the  motion  of  the  vessel. 

Thirty-one  data  records,  each  about  twenty  to  thirty  minutes  long,  were  recorded 
during  the  test  period  of  phase  one.  Eight  of  these  data  sets  were  reduced  to  power 
spectral  densities  (PSD)  in  terms  of  the  wave  encounter  frequencies.  Each  PSD  value 
has  about  fifty  degrees  of  freedom. 

The  results  show  the  significant  wave  height  (SWH)  as  measured  by  the  two  radars 
to  agree  within  about  ten  percent.  The  visual  estimates  of  the  SWH  show  a  larger 

a8reement/elativc  to  the  radar  derived  SWH.  One  buoy  measurement  was 
available  to  compare  with  the  radar  results.  The  buoy  derived  SWH  (0.5m)  was  about 
one-half  the  radar  measurement.  This  is  due  to  a  light  off  shore  wind  that  day  that 

forUJh«daln  f“n  SWH.for  both  the  bu°y  *nd  the  radars  as  well  as  a  shorter  fetch 
for  the  development  of  waves  at  the  buoy  site  as  conmarerf  w-it-h  t-ka  .1.4.  1 _ _ _ v 
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Introduction: 

(References  in  the  following  text  are  listed  at  the  end  of  FORWARD) 

The  U.  S.  Coast  Guard  (KSCG)  is  interested  in  continuing  research 
on  the  springing  phenomena  experienced  by  Great  Lakes  ore  carriers  and  on 
the  effect  springing  stresses  have  on  the  hull  girder  longitudinal 
strength.  A  vessel's  structure  is  very  complex  when  it  comes  to 
analyzing  the  structural  response  to  static  (cargo,  fuel,  ballast,  etc.) 
and  dynamic  (slarming,  wave-induced,  springing)  loads.  Of  primary 
interest  is  the  combined  effect  of  the  springing  and  wave  induced 
stresses  (Figure  I). 


FIGURE  I 

Springing  is  induced  by  encountered  waves  matching  the  natural 
frequency  of  a  vessel  which  creates  a  resonant,  vertical  2-node 
oscillation  of  the  hull  girder.  At  times,  springing  may  develop  1/3  to 
1/2  of  the  total  stress  level  within  the  hull  girder.  This  is  one  area 
yet  to  be  completely  verified.  The  longest  ships  plying  the  Lakes  have 
reached  1000  feet  with  an  L/D  ratio  of  20/1  to  22/1.  These  vessels  are 
very  long  and  flexible  and  are  susceptible  to  the  effects  of  combined 
wave-induced  and  springing  stresses.  One  such  vessel  is  the  M/V  STEWART 
J.  CORT. 
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Mr.  Don  Hammond  of  the  Advanced  Space  Sensing  Application  Branch, 
Naval  Research  Laboratory  (NRL),  has  submitted  this  Final  Report,  “Great 
Lakes  Wave  Height  Radar  System",  as  the  preliminary  evaluation  of  the 
wave  measurement  system  installed  on  the  M/V  STEWART  J.  CORT.  Of  primary 
importance  during  the  '78  and  '79  operating  seasons  of  the  CORT  has  been 
to  test  and  evaluate  the  wave  measurement  system.  The  discussion  below 
will  provide  a  description  of  how  this  research  fits  into  the  overall 
USCG  Great  Lakes  springing  program. 

Past  Research: 

The  CORT  has  been  the  subject  of  past  research  with  the  full 
cooperation  and  support  of  her  owner  and  operator,  Bethlehem  Steel 
Corporation.  Teledyne  Engineering  Services  (TES)  (formerly:  Teledyne 
Materials  Research  Corporation),  under  contract  to  the  USCG,  instrumented 
the  CORT  between  1971  and  1974  to  study  the  bending  stresses  experienced 
by  the  vessel  during  normal  operations.  The  vessel  was  instrumented  with 
five  sets  of  strain  gauge  transducers,  a  metritape  liquid-level  sensor  to 
measure  wave  encounter  periods,  bow  and  stern  vertical  accelerometers,  a 
pendulum  located  near  the  center  of  gravity  for  measuring  pitch,  and  a 
data  recording  system.  More  complete  details  on  this  research  and 
Teledyne's  analysis  can  be  found  in  Reference  A. 

During  the  period  of  1973  to  1975,  Mr.  Dalzell  of  Davidson 
Laboratory,  Stevens  Institute  of  Technology,  was  contracted  by  the  Ship 
Structure  Committee  (SSC)  to  reduce  and  analyze  the  wave  meter  data 
obtained  on  the  Sea-Land  McLean  and  to  perform  a  comparison  of  the  two 
wave  measurement  devices  used  on  the  vessel;  specifically,  a  microwave 
radar  and  a  Tucker  Meter.  Mr.  Dalzell  concluded  that  "The  evidence 
strongly  suggests  that  neither  of  the  wave  measuring  systems  can  be 
regarded  as  a  standard  by  which  the  other  may  be  judged."  (Ref.  B,  Vol. 

10)  Reference  B  describes  this  work  in  more  detail. 

USCG  Springing  Research  Program: 


In  an  effort  to  build  upon  this  past  research,  the  USCG  has 
undertaken  a  four  part  program  to  obtain  a  more  thorough  understanding  of 
springing. 

Part  1.  Wave  Measurement  Device  Survey  -  (Reference  C) 

In  September  of  1976,  NRl  was  contracted  by  the  USCG  to  survey 
existing  wave  measurement  devices  capable  of  providing  wave  height  and 
direction.  NRL  took  into  consideration  the  conclusions  and 
recommendations  reached  by  Mr.  Dalzell  (Ref.  B).  As  a  result,  NRL 
recommended  (1)  to  replace  the  marine  radar  on  the  CORT  with  a  Raytheon 
Model  1020  marine  radar  which  would  provide  the  wave  images  necessary  for 
determination  of  wave  length  and  wave  direction  via  a  laser  fourier 
transform  using  an  optical  processor  and  (2)  to  install  a  Collins  Radio 
Co.  radar  altimeter  for  measuring  wave  height.  The  USCG  considered  these 
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recommendations  and  finally  authorized  a  modified  wave  measurement  system 
to  be  installed  on  board  the  CORT.  This  system  is  described  in  Part  3 
below  and  within  the  text  of  this  report  by  NRL. 

Part  2.  Damping  Coefficient  Determination  -  (Reference  D) 

During  a  three  year  period  (1972-1975),  stress  measurements  were 
collected  from  three  Great  Lakes  ore  carriers,  namely  the  M/V  CORT,  the 
M/V  ROGER  BLOUGH,  and  the  S/S  RYERSON.  (References  F,G  and  H)  In  June 
1R77,  TES  was  contracted  by  the  USCG  to  transfer  specified  analog  stress 
records  from  these  3  vessels  to  FM  magnetic  tape.  These  stress  records 
contained  intervals  representing  high,  medium,  and  low  stress  levels 
during  loaded  and  ballasted  conditions.  These  magnetic  tapes  and  the 
associated  quick-looks  were  delivered  to  Structural  Dynamic  Research 
Corporation  (SDRC),  under  separate  USCG  contract,  to  determine  the  system 
damping  (combined  hydrodynamic,  structural  and  cargo  damping)  coefficient 
for  each  specified  interval  of  the  stress  records.  An  evaluation  of 
these  coefficients  was  performed  by  the  USCG  to  determine  a  relative 
accuracy  of  the  damping  coefficients  determined  by  SDRC. 

SDRC  used  an  analytical  technique  known  as  MODAMS  (MODal  Analysis 
and  Modeling  System).  The  technique  used  to  determine  the  modal 
parameters  desired  was  the  multi-degree-of-freedom  (MDOF)  curve  fit  to 
the  frequency  response  data  over  a  specified  frequency  range.  The  USCG 
evaluation  compared  SDRC's  coefficients  with  coefficients  derived  f'-om  a 
mechanical  analysis  of  independent  stress  records  from  these  vessels. 

The  damping  coefficients,  the  damping  frequency  and  the  logarithmic 
decrement  were  used  in  the  comparison.  The  evaluation  concluded  that 
SDRC's  results  were  reasonable  but  not  totally  conclusive.  Further 
damping  definition  may  be  provided  via  full  scale  hull  exciters.  Further 
detail  on  SDRC's  analysis  and  the  evaluation  may  be  found  in  Reference  D. 

Part  3.  Wave  Measuring  Device  T&E  -  (This  Report) 

In  December  1977,  WL  was  contracted  by  the  USCG  to  install  a  wave 
measurement  system  and  to  perform  a  test  and  evaluation  (T&E)  of  this 
system.  The  wave  measurement  system  proposed  by  NRL  (Part  1)  was 
modified  to  measure  wave  height  from  two  different  radars  and  to 
innediately  input  the  information  via  an  analog/digital  interface  to  an 
on-board  data  acquisition  system.  The  decision  was  made  to  opt  on  the 
side  of  obtaining  a  point  spectra  (in  lieu  of  a  directional  spectra)  with 
a  PDP-11/03  data  acquisition  system  rather  than  go  for  a  directional 
spectra  with  no  acquisition  system  capable  of  on  board  data  analysis. 

The  wave  measurement  system  consists  of  a  NRL  microwave  radar  and  a 
Collins  Radio  Co.  radar  altimeter,  two  sets  of  accelerometers  for 
reduction  of  ship  motions,  and  a  roll  potentiometer. 

Validation  of  the  on-board  measurements  was  attempted  by  comparing 
these  with  wave  measurements  of  a  NOAA  wave  buoy  positioned  2  miles 
offshore  of  Copper  Harbor,  Michigan  in  Lake  Superior.  This  effort  proved 
to  be  unsuccessful  due  to  (1)  infrequency  of  rendezvous  with  the  buoy, 

(2)  distance  of  CPA  with  buoy  (varied  from  2  to  10  miles)  and  (3)  the 
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difference  In  water  depth  and  wind  fetch  between  the  buoy  and  the  CORT. 
Wl's  report  evaluates  the  data  obtained  by  comparing  spectral  analyses 
of  both  radars  with  visual  observations  and  with  a  few  NOAA  buoy 
measurements.  Generally,  the  comparisons  show  some  agreement  but  the 
evidence  is  not  conclusive  and  stills  raises  many  questions  and  many 
doubts. 

Part  4(a)  Full  Seaway  Instrumentation  -  (In  progress  during  '79  and  ’80) 

This  project  is  an  extension  of  Part  3  to  measure  hull  girder 
stresses  concurrently  with  measuring  desired  wave  parameters  on-board  the 
M/v  STEWART  J.  CORT  during  normal  operations  in  fiscal  year  1979.  The 
results  of  this  data  collection  project  will  provide  stress  and  wave  data 
measured  in  a  simultaneous  manner  that  can  be  used  by  later  investigators 
to  explain  the  phenomenon  of  springing  and  wave  induced  responses  as  well 
as  the  method  for  combining  the  two.  With  the  on-board  data  acquisition 
and  reduction  system  on  the  CORT,  many  of  the  questions  related  to  these 
phenomena  may  be  investigated  during  this  full  scale  instrumentation 
project.  This  work  will  be  performed  by  David  Taylor  Naval  Ship  Research 
and  Development  Center  (DTNSRDC). 

Part  4(b).  Full  Scale  Pressure  Distribution  Measurements  -  (In  progress 
during  '79  and  '80) 

As  an  add-on  to  Part  4,  the  Ship  Structure  Comnittee  (SSC)  has 
contracted  with  DTNSRDC  to  collect  full  scale  pressure  distribution 
measurements  concurrently  with  the  wave  and  stress  measurements  on  the 
M/v  CORT.  Fifteen  50  psig  pressure  transducers  have  been  installed  in 
the  forward  section  of  the  vessel.  These  pressures  are  scheduled  to  be 
analyzed  by  the  American  Bureau  of  Shipping  (ABS)  in  conjunction  with  a 
time-domain  analysis  of  the  wave  heights. 

Figure  II  shows  the  instrumentation  on  the  CORT  that  will  be  used 
for  full  scale  measurements  during  the  '79  -  '80  operating  season. 
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The  validation  of  the  wave  measurement  system  remains  foremost  in 
importance  as  the  full  scale  measurements  continue.  During  the  1979  fall 
season,  DTNSRDC  will  continue  the  full  scale  measurements  on  the  CORT. 
This  time  validation  of  the  wave  measurement  system  will  be  attempted  by 
deploying  a  wave  buoy  approximately  2  miles  ahead  of  the  vessel  via  USCG 
helicopters.  Helicopters  were  available  from  USCG  Air  Station  Chicaqc 
(Illinois)  and  USCG  Air  Station  Traverse  City  (Michigan).  Two  modified 
Eastech  Model  226  wave  buoys  will  be  utilized  for  this  project. 

i 

Bethlehem  Steel  Corporation,  the  Master  of  the  CORT,  the  ships'  crew, 
and  the  personnel  at  USCGAS  Chicago  and  USCGAS  Traverse  City  have 
provided  their  full  cooperation  and  support  for  this  research  project. 


Mark  D.  Noll,  Project  Officer 
LT,  U.S.  Coast  Guard 
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INTRODUCTION 


The  ore  vessel  M/V  STEWART  J.  CORT  was  the  first  of  the 
Great  Lakes  vessels  with  a  length  as  great  as  one  thousand 
feet.  For  this  reason,  it  has  been  a  subject  of  particular 
interest  for  the  study  of  hull  vibration  and  springing. 

One  of  the  early  studies  (reference  1)  used  five 
vertical  accelerometers,  one  on  each  side  near  the  bow,  one 
at  the  quarter  point,  one  at  the  mid  point  and  one  at  the 
stern,  and  a  midship  strain  gage.  Based  on  the  analysis  of 
records  taken  in  Lake  Superior  when  the  vessel,  inQthe 
loaded  condition,  encountered  4  to  6  foot  3eas,  60  off  the 
bow,  a  determination  was  made  of  the  maximum  first  mode 
deflection  profile,  the  RMS  deflection  profiles  for  the 
first  three  modes  of  flexural  vibration,  and  the  maximum  and 
RMS  midship  stress  values.  An  evaluation  of  the  first  mode 
damping  was  unsuccessful  because  the  Cort  was  undergoing 
springing  in  the  stationary  anchor  drop  position. 

Another  study  (reference  2)  analyzed  the  midship  bending 
stress  obtained  from  a  deck  mounted  strain  gage  for  two 
half-hour  intervals  of  high  stress  for  the  M/V  STEWART  J. 

CORT  and  the  SS  CHARLES  M.  BEEGHLY.  The  stress  records  were 
analyzed  for  the  RMS  springing  and  wave  bending,  the  maximum 
peak  values,  and  the  instantaneous  values. 

The  present  study  is  the  first  phase  of  a  two-phase 
in-depth  study  of  vessel  springing  induced  by  surface  waves. 
This  first  phase  involves  the  installation  of  a  wave  measure¬ 
ment  system  on  the  M/V  STEWART  J.  CORT,  the  recording  of 
wave  data  and  the  evaluation  of  the  wave  measurement  system. 
The  second  phase  will  be  done  during  the  year  1979  and  will 
include  the  simultaneous  recording  of  the  waves  encountered 
by  the  CORT  and  the  response  of  the  vessel  as  recorded  by 
hull-mounted  strain  gages. 
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WAVE  MEASUREMENT  SYSTEM 


The  wave  measurement  system  consists  of  two  radars 
pointed  toward  the  water  surface  forward  of  the  vessel.  Each 
radar  has  an  associated  pair  of  accelerometers  to  record  the 
motion  of  the  vessel.  The  two  radar  systems  were  available 
from  earlier  programs  of  the  Ship  Structure  Committee  and 
the  U.S.  Coast  Guard.  One  an  NRL  radar,  was  built  and  used 
on  the  studies  of  the  SL-7  class  container  ship.  References 
(3,  4,  5)  describe  the  former  SL-7  program.  The  other,  a 
Collins  radar,  was  used  to  assess  its  value  in  predicting 
slamming.  The  use  of  two  radars  serves  to  provide  a  compari¬ 
son  of  results  and  also  provides  redundant  reliability. 


The  system  configuration  was  slightly  changed  after  the 
end  of  the  first  phase  of  operations.  The  Collins  radar 
horns,  along  with  two  new  motion  accelerometers,  were  mounted 
on  a  new  boom  over  the  bow  to  reduce  any  effect  the  dynamic 
bow  wake  might  have  on  the  Collins  radar  return.  The  stable 
platform  that  held  a  pair  of  accelerometers  and  a  roll  sensor 
were  removed  from  the  pilot  house  top.  Repeated  efforts 
throughout  phase  one  were  unsuccessful  in  getting  the  stable 
platform  to  work.  During  phase  two  a  roll  gyro  will  be  used 
to  record  the  roll  angle  of  the  vessel. 

Figure  1  shows  the  locations  of  the  two  radar  antennas 
on  the  vessel  and  their  footprints  on  the  water  surface  in 
front  of  the  vessel's  bow.  The  NRL  radar  antenna  is  located 
on  the  forward  port  corner  atop  the  pilot  house.  Its  3  1/2 
degree  wide  antenna  beam  is  directed  forward  at  a  vertical 
angle  of  25  degrees.  The  footprint  illuminated  on  the  sur¬ 
face  is  nearly  circular  with  1.2  meter  diameter.  The  Collins 
radar  antennas  are  mounted  at  the  end  of  a  4.5  meter  long 
boom  directly  forward  of  the  bow  at  the  main  deck.  The  20 
degree  wide  beam  of  the  Collins  radar  antennas  are  directed 
forward  at  a  25  degree  vertical  angle.  The  Collins  radar 
footprint  has  a  4.75  diameter  footprint.  Because  of  the 
footprint  locations,  the  Collins  radar  records  waves  from 
all  points  around  the  bow,  whereas  the  NRL  radar  records 
waves  mainly  from  points  off  the  port  bow. 


The  NRL  radar  achieves  its  range  resolution  with  a  two 
nanosecond  (10  )  wide  pulse.  A  Tektronix  sampling  scope  is 

used  to  transform  the  narrow  video  pulses  to  equivalent  time 
pulses  that  can  be  processed  with  standard  logic  circuits. 
The  range  tracker  circuit  updates  pulse  to  pulse  the  pro¬ 
grammable  time  delay  that  is  a  feature  of  the  scope.  The 
ra2ge^racker  senses  changes  in  return  pulse  time  position 

in deiay  t0  keeP  th*  video  pulse  Centered 
in  a  15  meter  range  window.  The  Collins  radar  uses  a 

°fFits  “.neuter  frequency  to  oht Jn  renge  " 

solution.  Figure  2  shows  the  sawtooth  shaped  modulation 
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RADAR  ANTENNA  LOCATIONS  ON  T|ir  VESSEL 


waveform.  Figure  2  also  shows  that  the  return  waveform  is 
shifted  in  time,  due  to  the  round  trip  time  of  flight  of  the 
signal.  The  Collins  radar  uses  the  instantaneous  frequency 
difference  to  measure  distance  in  this  fashion. 

The  NRL  radar  signal  is  in  digital  binary  coded  format 
for  direct  recording  by  the  PDP-11/03  computer.  The  samples 
of  NRL  radar  distance  are  sent  to  the  computer  ten  times  a 
second.  The  raw  counts  are  converted  to  meters  by  multi¬ 
plication  by  0.01875.  This  factor  is  derived  from  the  fact 
that  each  count  represents  exactly  1/8  of  a  nanosecond  round 
trip  delay.  Each  nanosecond  in  turn  is  equivalent  to  0.15 
meters  distance.  Therefore,  each  count  is  1/8  times  0.15 
meters  or  0.01875  meters. 

The  Collins  radar  was  calibrated  by  noting  the  change 
in  the  digitized  counts  recorded  when  a  known  length  of  co¬ 
ax  cable  was  added  between  the  receiver  input  and  the  receiver 
horn  antenna.  The  length  of  the  additional  transmission 
line  was  8.5344  meters.  The  increase  in  the  average  count 
was  524.44.  The  ratio  of  the  speed  of  the  radar  signal  in 
air  to  that  in  the  cable  is  1.492:1.  This  gives  a  calibration 
factor  of  0.01214  to  convert  Collins  radar  counts  to  meters 
distance. 

The  appendices  C,  D,  and  E  contain  descriptions  of  both 
the  NRL  and  Collins  radar  systems. 


SYSTEM  CONFIGURATION 

The  signals  from  the  sensing  devices  were  all  brought 
by  cable  to  the  signal  distribution  chassis  and  then  to  the 
PDP-11/03  computer  during  phase  one  as  shown  in  figure  3. 

The  NRL  radar  signal  is  in  twelve  bit  digitial  format.  The 
NRL  radar  provides  a  data  ready  flag  ten  times  a  second  to 
serve  as  an  interrupt  to  the  computer  program.  The  computer 
program  stores  the  NRL  digital  range  signal  and  then  it  uses 
its  analog-to-digital  converter  to  sample  the  analog  signals 
from  the  other  devices.  The  program  continues  in  this 
manner  for  the  time  duration  entered  by  the  operator  in  the 
header  log.  The  Collins  radar  signal  strength  channel 
changes  to  a  large  value  if  the  Collins  return  signal 
strength  is  below  the  valid  signal  threshold. 

The  header  log  is  produced  by  the  computer  operator 
responding  via  the  Decwriter  to  each  parameter  listed  by  the 
computer.  the  computer  program  will  cause  the  heading,  date 
and  time  to  be  typed  followed  by  a  standard  list.  Table  1 
shows  a  typical  phase  one  header  log.  There  is  a  provision 
at  the  bottom  for  the  operator  to  type  in  three  comment 
statements,  each  limited  to  sixty  characters. 
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FREQUENCY 


COLLINS  RADAR  MODULATION  WAVEFORM 


FIGURE  2 
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FIGURE  3 
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TABLE  I 


SERIAL  NO.  29 


GREAT  LAKES  WAVE  HEIGHT  PROJECT 
DATE:  16  NOV  78  TIME: 

DURATION  OF  RUN  IN  MINUTES  IS 
NORTH  LATITUDE  (DD  MM) 

WEST  LONGITUDE  (DD  MM) 

VESSEL'S  SPEED  (MPH  -  XX. X) 

VESSEL'S  HEADING  (DEGREES) 

VESSEL'S  DRAFT  (FEET) 

WIND  DIRECTION  (DEGREES) 

WIND  SPEED  (KNOTS) 

WAVE  DIRECTION  (DEGREES) 

WAVE  HEIGHT  (FEET) 


22:52:33 

4 

47  29 
87  54 
16.3 
270 
17 
315 
20 
315 
4 


REMARK  LINE  -  UP  TO  60  CHARACTERS 
COPPER  HARBOR  1.5  MILES,  CODE  3  F  17,  M  18-11,  AFT  21-3 


REMARK  LINE  -  UP  TO  60  CHARACTERS 
JUST  WEST  OF  BUOY 


REMARK  LINE  -  UP  TO  60  CHARACTERS 


CHANGED  HEAD  TO  265 


Figure  4  is  a  flow  diagram  of  the  computer  program  to 
record  data.  The  header  is  followed  by  the  recording  of  the 
NRL  radar  sample  taken  at  the  data  ready  flag  time.  This  is 
followed  by  a  sequential  scan  by  the  computer  analoq-to- 
digital  converter  of  the  eight  analog  signal  channels.  The 
eighth  analog  channel  is  not  connected  for  this  phase  of  the 
project  and  is  regarded  as  a  spare. 


ACCELEROMETER  CALIBRATION 

The  accelerometers  are  mounted  in  pairs  to  sense  the 
vertical  movement  of  the  vessel.  The  accelerometers  were 
calibrated  by  the  usual  method  of  recording  the  static 
output  voltage  as  the  angular  position  of  the  accelerometer 
is  varied  relative  to  the  gravity  vector.  Figure. 5  is  the 
plot  of  the  calibration  data  for  the  unstabilized  vertical 
accelerometer  mounted  on  the  NRL  radar  mount  on  top  the 
pilot  house.  A  least  squares  fit  to  the  calibration  data 
gives  a  factor  of  9122  counts  per  "g".  Figure  6  is  the  plot 
of  the  unstabilized  horizontal  accelerometer  calibration 
data. 


The  Setra  System,  Inc.,  Model  100  accelerometers  used 
on  the  CORT  have  a  high  natural  frequency  of  350  hertz. 

This  is  so  far  above  any  frequency  at  which  the  vessel  can 
spring  that  the  static  calibrations  are  valid  for  all  fre¬ 
quencies  observed. 

Two  Setra  System,  Inc. ,  Model  114  accelerometers  were 
mounted  on  the  antenna  end  of  the  new  longer  bow  boom  for 
the  Collins  Radar.  The  calibration  factor  for  the  vertical 
boom  accelerometer  is  5325  counts  per  "g". 

The  sense  of  the  vertical  accelerometer  is  such  that  a 
downward  acceleration  results  in  a  reduction  in  the  counts 
recorded  and  an  upward  acceleration  is  recorded  as  an  increase 
in  the  recorded  counts. 

Thus,  there  are  two  sets  of  accelerometers  for  use 
during  the  phase  two  operations.  One  set  is  mounted  at  the 
NRL  radar  mount  and  the  other  set  is  mounted  at  the  Collins 
radar  antenna  boom  over  the  vessel's  bow. 


SYSTEM  CALIBRATION  CONSTANTS 

The  conversion  of  the  raw  counts  recorded  from  each 
device  in  the  wave  measurement  system  to  engineering  units 
is  done  by  multiplication  by  the  following  factors: 
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COMPUTER  FLOW  CHART 
DATA  RECORDING  PROGRAM 

FIGURE  4 
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0  VERTICAL  ANGLE 


CALIBRATION  VERTICAL  UNSTABILIZED  ACCELEROMETER 

17  NOV  1978 


FIGURE  5 


CALIBRATION  HORIZONTAL  UNSTABILIZE 0  ACCELEROMETER 

17  NOV  1978 

FIGURE  6 
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DEVICE  FACTOR  ENGINEERING  UNIT 

PILOT  HOUSE  TOP 


NRL  Radar  Range 

0. 01875 

Me  te  r  s 

Vertical  Accelerometer 

1.096  x 

10"4 

g's 

Horizontal  Accelerometer 

1.938  x 

io-4 

g’ s 

END  OF  BOW  BOOM 

Collins  Radar  Range 

0.01214 

Meters 

Vertical  Accelerometer  * 

1.783  x 

10"4^ 

g's 

Horizontal  Accelerometer  * 

2.520  x 

10-4 

g  '  s 

*  For  Phase  two  measurements  only 


SPATIAL  FILTER  EFFECT  OF  THE  RADAR  FOOTPRINT 

Since  the  two  radars  measure  the  distance  to  an  area  on 
the  water  surface  it  is  necessary  to  analyze  how  the  diameter 
of  the  radar  footprint  affects  the  results.  It  is  reasonable 
to  assume  that  the  amplitude  of  a  water  wave  with  a  crest- 
to-crest  distance  that  is  short  compared  with  the  diameter 
of  the  radar  footprint  will  be  "averaged  out".  That  is,  the 
radar  distance  will  be  the  result  of  a  spatial  average  over 
the  footprint.  One  dimensional  analysis  of  this  effect 
is  adequate  because  the  radar  illumination  function  is 
symetrical  about  its  center  axis. 

The  illumination  function  for  each  radar  due  to  the 
directivity  of  its  antenna  is  a  gaussian  function  of  angle 
relative  to  its  center  axis. 

fi2 

D  ( 0)  =  A  exp-— 5-  (1) 

B 

"A"  is  a  factor  such  that  the  definite  integral  of  D  is  one. 
Where  0  is  the  angular  direction  relative  to  the  beam  center 
axis  and  B  is  a  measure  of  the  width  or  fall  off  of  intensity 
with  angle.  For  small  angles,  0  is  approximately 

8  ■  I  '2> 

Where  0  is  measured  in  radians,  S  is  the  distance  on  the 
surface  from  the  point  of  intersection  of  the  surface  with  the 
antenna  beam's  center  axis  and  H  is  slant  range  of  the 
antenna  to  the  surface. 
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This  substitution  for  6  in  equation  (l)  gives 


D  (S )  =  ~~  exp  '  (3) 

JttHB  H  B 

A  water  wave  component  with  wave  number  k  will  have  its 
apparent  amplitude  reduced  as  a  result  of  being  averaged 
with  equation  (3)  as 

(+oo 

D  (S)  cos  (kS+'l' )  dS  (4) 


Where  k  is  2tt  divided  by  the  wavelength  of  the  water  wave  X 
and  Y  is  the  phase  of  the  wave  component  relative  to  the 
radar  beam  center  axis. 

The  result  of  scanning  the  cosine  shaped  water  wave  of 
wave  number  k  is  to  reduce  its  amplitude  by  the  following 
factor. 


F  (k) 


2  2  2 
,k  HZB  , 
exp-  ( - 2 - ) 


(5) 


The  relation  used  to  find  the  critical  wave  number  (k  ) , 
that  is  where  the  variance  is  reduced  by  one  half  dueCto 
spatial  filtering,  is  found  by  setting  the  exponent  function 
equal  to  the  natural  logarithm  of  the  square  root  of  two  and 
solving  for  kc-  This  can  then  be  related  to  the  encounter 

frequency  for  different  encounter  angles. 


(6) 


(7) 

fe  =  >fgkc  +  Ukc  cost))  (8) 

C  ~Ti t  SW 


2  2  2 

k  Vb  _  1  - 

c  —  -re  ion  2 
- 4  2 


v  =  (2  £n  2) 
C  HB 


1/2 


Where  fe  is  the  frequency  of  encounter  in  Hertz,  g  is 
the  acceleration  of  gravity,  U  is  the  speed  of  the  vessel, 
and  <J>  is  the  angle  between  the  vessel’s  heading  and  the 
water  waves.  The  speed  of  the  vessel  is  very  close  to  6.7 
meters/second  (15  miles  per  hour)  for  all  data  sets.  Next 
this  relation  is  applied  to  both  radars  as  follows. 
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The  illumination  function  for  the  Collins  radar  is 
derived  by  the  relation  between  the  half  power  angular  width 
or  field  of  view  (FOV)  and  the  antenna  gain  factor.  The 
Collins  radar  antenna  horns  have  a  gain  of  100.  Since 

FOV  =  arc  cos  (l-£) ,  (9) 

yj 


the  angular  FOV  is  0.2  radian  for  the  one  way  antenna  beam. 
The  two  way  beam  width  is  narrower  by  the  factor  of  the 
square  root  of  one  half. 

From  ( 1 ) : 


D ( . 1414)  =  j  =  A  exp 
en2  . 


B 


,VI 


B 


10 (£n2) 


172 


(10) 

(ID 

(12) 


B  =  0.16986  radians 


The  slant  range  (H)  of  the  Collins  radar  antenna  from  the 
surface  is  12  meters. 


The  NRL  radar  has  a  60  centimeter  diameter  parabolic 
antenna.  The  NRL  radar  wavelength  is  3  centimeters.  This 
combination  gives  a  one  way  antenna  angle  of  0.0611  radians. 
The  two  way  antenna  half  power  beam  width  is  0.0611A/7. 
Therefore,  the  beam  factor  B  for  the  NRL  radar  is 

B  =  (.0611) 

^2  in  2 

B  =  .05189  radians 


The  slant  range  of  the  NRL  radar  antenna  (H)  from  the  surface 
is  23  meters. 

Table  2  is  a  list  of  the  encounter  frequencies  (fe)  for 
which  the  PSD  is  reduced  to  one  half  its  true  value  due  to 
the  spatial  averaging  of  the  radar's  illumination  function. 
The  angle  <J>  is  zero  for  seas  head  on  the  bow.  The  angle 
increases  as  the  seas  are  off  either  the  starboard  or  port 
bow.  For  seas  close  on  the  bow  the  Collins  spatial  filter 
begins  to  reduce  the  variance  of  the  encounter  frequency  for 
components  above  one  Hertz.  The  NRL  radar  spatial  filtering 
becomes  a  significant  factor  for  encounter  frequencies  above 
1.5  Hertz.  The  data  shows  the  peak  in  the  PSD  to  fall 
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TABLE  2 


<J>  Degrees 

Encounter  Frequency  for 

1/2  Variance  Response 

OFF  THE  BOW 

COLLINS 

NRL 

HERTZ 

HERTZ 

0 

0.994 

1.547 

10 

0.985 

1.531 

20 

0.957 

1.483 

30 

0.912 

1.406 

40 

0.851 

1.300 

50 

0.775 

1.171 

60 

0.687 

1.020 

70 

0.589 

0.855 

80 

0.486 

0.678 

90 

0.379 

0.495 

between  0.3  and  0.4  Hertz.  Therefore,  spatial  filtering  by 
either  radar's  antenna  pattern  should  not  affect  the  results. 
This  applies  to  head  to  beam  seas  only.  Following  seas  are 
more  complicated  and  are  not  addressed  because  the  main 
interest  of  this  study  is  head  seas. 


BOW  WAKE 

The  bow  of  the  M/V  STEWART  J.  CORT  is  a  blunt,  almost 
vertical  wall  with  a  small  amount  of  slant  to  either  side. 

The  fifteen  miles-per-hour  speed  of  the  vessel  causes  the 
bow  to  exert  a  great  deal  of  force  on  the  water  in  front  of 
it.  This  additional  pressure  raises  the  level  of  the  water 
just  in  front  of  the  vessel  above  the  average  surface  of  the 
lake.  This  forward  bulge  of  water  extends  perhaps  to  one- 
fifth  the  length  of  the  CORT  or  about  60  meters  forward. 

This  standing  wave,  relative  to  the  vessel .appeared  to  have 
very  little  influence  on  the  wave  measurements  for  the  condi¬ 
tions  encountered  during  these  field  trials. 

There  is  a  dynamic  or  time  varying  aspect  to  the  bow 
wake  that  could  show  up  in  the  PSD  at  very  low  frequencies. 
This  dynamic  portion  of  the  bow  wake  is  a  surge  of  water 
that  seems  to  explode  periodically  just  in  front  of  the 
vessel.  It  looks  similar  to  the  surf  on  a  beach.  The 
pressure  appears  to  build  just  below  the  surface  and  is 
released  as  an  upward  and  forward  surge  of  water  with  much 
foam.  It  runs  down  and  is  quiet  for  a  time  and  repeats 
again.  The  size  and  timing  depend  on  the  relative  direction 
and  size  of  the  water  waves.  It  is  larger  for  head  seas  and 
large  wave  heights.  One  twenty-minute  period  of  visual 
observations  with  head  seas  and  moderate,  1.5  meters,  wave 
heights  showed  that  the  foam  extended  as  far  forward  as  7.5 
meters. 

These  visual  observations  led  to  the  installation  of  a 
4.5  meter  long  boom  to  hold  the  Collins  radar  horns  out  in 
front  of  the  dynamic  bow  wake  for  phase  two  measurements.  The 
Collins  radar  horns  are  tilted  forward  with  a  vertical  angle 
of  25  degrees,  thus  placing  their  footprint  9.7  meters  forward 
of  the  bow.  Figure  7  shows  the  relative  size  and  location  of 
the  Collins  radar  for  phase  one  and  phase  two  operations. 

The  new  boom  is  a  triangular  truss  structure,  30  centi¬ 
meters  on  a  side.  It  is  a  standard  structure  for  a  radio 
antenna  tower.  It  is  a  very  stiff  structure  made  of  steel. 

A  guy  wire  near  the  front  end  adds  additional  support.  It  is 
mounted  on  the  top  of  a  steel  pipe  standing  on  the  main  deck 
at  the  bow.  It  is  rotated  in  on  the  port  side  for  storage 
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COLLINS  RADAR  FIELD  OF  VIEW 


between  runs  or  while  going  through  the  Soo  locks.  The 
motion  sensing  accelerometers  are  mounted  in  a  weather  proof 
can  right  at  the  forward  end  of  the  boom.  They,  therefore, 
measure  the  motion  of  the  Collins  horns  directly. 


DATA  RECORDED 

There  were  31  floppy  disks  of  data  recorded  during  the 
October  November  1978  period.  Table  3  is  a  serial  list  of 
the  data  by  date  with  a  notation  of  the  significant  wave 
heights.  The  data  that  has  its  Power  Spectral  Density  plotted 
in  this  report  is  noted  by  an  asterisk.  The  header  log 
sheets  for  these  sets  are  included  in  Appendix  B.  The 
selection  of  the  data  was  based  on  (1)  close  approach  to 
Copper  Harbor  NOAA  wave  buoy,  (2)  Head  seas,  (3)  vessel 
springing,  and  (4)  large  wave  heights. 

Serial  numbered  runs  6,  7,  23,  27  and  29  were  near 

Copper  Harbor.  Serial  numered  runs  10  and  14  were  nearly 

head  seas.  Serial  numbered  runs  6,  7,  10,  14,  and  30  show 

springing  of  the  vessel  as  evident  by  the  vertical  acceler¬ 
ometer  PSD  coefficient  near  the  vessel's  resonant  frequency 
of  0.32  Hertz.  Runs  14  and  30  have  large  wave  heights. 

It  is  important  to  note  that  the  accelerometer  PSD  of 
vertical  displacement  is  multiplied  by  100  so  that  it  can  be 
plotted  on  the  same  scale  as  the  radar  PSD  of  range.  In 
general,  the  accelerometer  shows  a  narrow  band  of  energy 
near  the  vessel's  springing  frequency  of  0.32  Hertz.  A 
comparison  of  the  maximum  estimated  peak-to-peak  deflections 
near  0.32  Hertz  of  run  serial  number  30  of  7.6  centimeters 
with  Critchf ield’ s  (reference  1)  maximum  of  11.1  centimeters 
shows  reasonable  agreement. 


DATA  ANALYSIS 

The  data  counts  are  reduced  and  analyzed  in  terms  of 
power  spectra  components  by  a  Fast  Fourier  Transform  (FFT) 
computer  program  written  by  Dave  Walden  of  the  U.S.  Coast 
Guard.  The  program  allows  the  operator  to  select  the  dura¬ 
tion  of  the  data  used  for  analysis.  The  numerical  printouts 
are  in  terms  of  total  sample  variance  at  a  particular  fre¬ 
quency  within  the  resolution  bandwidth. 

These  FFT  values  were  converted  to  engineering  units  of 
centimeters  squared  via  the  calibration  constants  and,  to 
make  all  the  data  comparable,  each  component  was  divided  by 
the  resolution  bandwidth  expressed  in  Hertz.  The  results 
are  plotted  versus  encounter  frequency  in  Hertz  (one  Hertz 
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TABLE  3 


DATA  RECORDS 


Serial 

Number 

Date 

1978 

Signif  i< 
Height-1 

COLLINS 

RADAR 

rant  Wave 
4eters 

NRL 

RADAR 

VISUAL 

BUOY 

1 

10-9 

2 

10-9 

3 

10-9 

4 

10-9 

5 

10-9 

6 

Ml)  (3) 

10-11 

1.5 

1.53 

7 

Ml)  (3) 

10-11 

1.36 

1.53 

8 

10-11 

9 

10-12 

10 

*  ( 2 )  (3) 

10-12 

2.5 

1.83 

11 

10-21 

12 

10-21 

13 

10-22 

14 

*  ( 2 )  (3) 

10-22 

2.7 

2.44 

15 

10-22 

16 

10-22 

17 

10-22 

18 

10-24 

19 

10-24 

20 

10-25 

21 

10-30 

22 

11-8 

23 

Ml) 

11-9 

1.7 

0.9 

24 

11-10 

25 

11-10 

26 

11-10 

27 

Ml) 

11-12 

1.1 

0.9 

0.52 

28 

11-13 

29 

Ml) 

11-16 

0.91 

0.8 

1.22 

30 

M3) 

11-23 

2.67 

2.98 

1.83 

31 

11-25 

1.96 

1.73 

1.53 

*  DATA  REDUCED 

(1)  Near  Copper  Harbor 

(2)  Head  Seas 

(3)  Vessel  Springing 
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equals  one  cycle  per  second) .  Therefore,  the  results  are 
plotted  in  the  form  of  Power  Spectral  Density  (PSD)  with 
engineering  units  of  centimeter  squared-seconds. 

Figure  8  is  a  simplified  flow  chart  for  the  FFT  program. 
The  operator  selects  one  of  the  data  files  (numbered  0  up  to  9) 
on  the  data  floppy  disk.  Next  the  number  of  data  samples  to 
be  used  for  the  FFT  analysis  is  selected.  The  next  option 
is  to  use  every  data  point  (ten  per  second)  or  every  other 
data  point  (five  per  second) .  The  operator  then  selects  the 
data  channel  for  analysis.  A  number  of  data  samnles  are 
selected  for  printout  in  the  event  the  operator  desires  to 
check  the  data.  The  program  prints  out  the  data  mean  value 
in  counts,  the  variance  and  root  mean  square  (rms)  values  in 
units  of  counts.  The  operator  then  has  the  option  of  using  a 
three  term  Blackman-Harris  window  function  on  the  time 
sampled  data  before  the  FFT.  Several  data  passes  were  made 
that  include  the  window  to  compare  with  results  without  the 
window.  The  window  was  not  used  on  any  of  the  data  plotted 
for  this  report  since  its  use  had  a  very  small  effect  on  the 
results.  The  next  step  is  to  indicate  how  many  FFT  values 
to  print.  The  operator  next  types  in  the  number  of  consecutive 
frequency  coefficients  to  sum.  Finally,  the  operator  types 
in  the  desired  number  of  coefficients  to  print  after  summing. 
These  printouts  are  the  total  variance  versus  frequency. 

The  data  analysis  is  based  on  procedures  for  the  analysis 
of  random  noise.  The  estimates  of  PSD  of  random  noise  have 
a  Chi-squared  probability  distribution.  The  stability  of 
the  estimates,  how  much  they  are  likely  to  change,  from  one 
data  segment  to  the  next,  depends  on  the  number  of  degrees 
of  freedom  in  each  estimate.  This  means  the  more  data  used 
to  obtain  a  single  PSD  coefficient  the  better  the  estimate. 

There  are  two  degrees  of  freedom  for  each  estimate  of 
the  PSD.  If  a  number  of  consecutive  values  are  summed  then 
the  degrees  of  freedom  are  increased  and  the  consistency  of 
the  estimate  is  improved.  This  improvement  is  offset  by  a 
poorer  frequency  resolution.  The  compromise  used  in  the 
plots  included  here  is  such  that  the  best  frequency  resolu¬ 
tion  possible  is  used  while  keeping  about  50  degrees  of 
freedom. 

The  header  logs  and  PSD  plots  for  8  sets  of  data  are 
included  in  this  report.  The  selected  data  includes  generally 
head  sea  conditions  or  close  approach  to  Copper  Harbor, 
Michigan.  NOAA  had  a  wave  buoy  near  Copper  Harbor,  Michigan 
and  some  of  that  data  is  used  for  comparison.  Some  of  the 
data  sets  show  the  vessel  springing  as  evident  in  the  accel¬ 
erometer  results. 
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SELECT 
INPUT  DATA 
FILE  AND 
CHANNEL 


SELECT 

NO.  OF  SAMPLES 
AND  SAMPLE  SPACING 


(See  Appendix  F 
for  detailed  fl 
diagram  and  cop 
of  program) 


FLOW  CHART 
DATA  REDUCTION 

FIGURE  8 
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The  plots  show  the  PSD's  of  the  Collins  Radar,  the  NRL 
radar  and  the  vertical  accelerometer  measurements  transformed 
to  displacement  by  double  integration.  Double  integration 
will  cause  the  derived  values  of  displacement  for  the  low 
frequencies  to  be  too  large.  Any  noise  in  the  acceleration 
measurement  will  be  emphasized  at  the  very  low  frequencies. 
The  unknown  constants  of  integration  also  show  up  as  very 
large  values  in  the  lowest  frequency  band.  Therefore,  the 
PSD  estimate  of  displacement  at  the  lowest  frequency  resolved 
should  be  ignored. 


VESSEL  MOTION  COMPENSATION 

A  large  ore  carrier  like  the  CORT  has  very  small  move¬ 
ments  in  pitch,  roll  and  heave  under  most  sea  conditions. 

The  corrections  to  the  radar  spectral  components  by  the 
accelerometer  measurements  are  likewise  smalx.  For  example, 
an  examination  of  the  radar  PSD  and  the  vertical  displace¬ 
ment  as  measured  by  the  double  integration  of  the  accelerom¬ 
eter  data  show  that  the  largest  displacement  PSD  is  about 
one  percent  of  the  radar  PSD.  The  amplitude  ratio  is  the 

square  root  of  the  PSD  ratio  or  ten  to  one. 

The  report  "Wavemeter  Data  Reduction  Method  and  Initial 
Data  from  the  SL-7  Containership”  (reference  5)  by  John 
Dalzell  uses  the  following  equation  to  relate  the  various 
components  of  encounter  spectrum. 

Sgg(We)  =  Srr(We)  +Szz(We)  +  2Crz(We) 

Sgg(We)  =  encounter  spectrum  of  waves 
Srr(We)  =  encounter  spectrum  of  corrected  radar  range 
Szz(We)  =  spectrum  of  heave  motion  via  the  accelerometer 
Crz(We)  =  Co-spectrum  of  radar  with  heave 

If  Szz(We)  is  very  small  compared  with  Srr(We)  then  Crz(We) 
will  also  be  small.  Thus,  the  approximation 

Sgg(We)  =  Srr(We) 

is  reasonably  good  for  this  case. 

The  PSD  plots  show  Szz (We)  to  be  about  one  percent  of 
Srr(We),  therefore,  both  Szz (We)  and  Crz(We)  are  small 
compared  with  Srr(We).  Therefore,  the  PSD  of  the  radar 
range  will  be  a  good  approximation  to  the  encounter  PSD  of 
the  vessel  and  the  waves. 
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DISCUSSION  OF  RESULTS 

A  comparison  of  the  significant  wave  heights  (SWH)  listed 
in  Table  3  show  the  visual  estimates  to  be  lower  than  the  radar 
estimate  five  times  and  higher  three  times.  The  small  SWH 
(less  than  about  one  meter)  have  the  largest  percent  difference 
in  the  estimates.  A  comparison  of  the  SWH  as  measured  by  the 
two  radars  for  the  three  cases  the  NRL  radar  was  working  show 
a  variation  of  about  ten  percent.  The  NRL  transmitter  was  not 
working  for  some  of  the  early  data  runs,  thus  of  the  eight 
sets  reduced  for  this  report,  only  runs  27,  29,  and  30  have 
both  radars  working. 

The  plots  of  Power  Spectral  Density  (PSD)  all  have 
about  50  degrees  of  freedom  for  each  estimate  which  says  that 
the  80%  confidence  interval  is  1.5  d£  wide.  They  are 
plotted  on  a  full  scale  of  10,000  cm  -second.  Larger  values 
of  full  scale  were  needed  to  display  numbers  14  and  30.  The 
values  of  PSD  of  the  vertical  displacement  calculated  from 
the  accelerometer  data  were  multiplied  by  100  so  that  they 
could  be  plotted  on  the  same  scale  as  the  radars.  The 
vertical  displacement  is  a  narrow  band  response  near  0.3 
Hertz  frequency.  The  amplitude  of  the  springing  in  runs  6, 

7,  10,  14,  23,  and  30  is  about  one  percent  of  the  variance 
of  the  PSD  of  the  waves  as  measured  by  the  radar  at  the  same 
frequency. 

Run  number  27  was  near  the  NOAA  wave  buoy  at  Copper 
Harbor,  Michigan.  The  buoy  data  was  transformed  into  en¬ 
counter  spectra  taking  into  account  the  observed  angle 
between  the  vessel  and  the  waves  as  recorded  on  the  log 
found  in  Appendix  B. 

Figure  9  is  a  map  of  the  region  around  Copper  Harbor. 

It  shows  the  vessel  heading  East  with  an  off  shore  wind  from 
the  South-East.  The  buoy  is  about  2  miles  from  the  shore 
and  the  vessel  is  about  5  miles  from  the  shore.  The  fact 
that  the  buoy  measures  a  half  meter  SWH  and  the  ship  radar 
indicates about  one  meter  SWH  is  due  to  the  longer  fetch 
(distance)  over  which  the  wind  blows  off  the  shore  to  the 
vessel.  Figure  10  is  a  plot  of  the  PSD  of  the  two  radars, 
the  buoy  data  transformed  to  encounter  spectra,  and  the 
vertical  displacement.  The  scale  has  been  reduced  to  better 
show  the  small  values  of  PSD. 

Table  4  is  the  NOAA  data  in  the  first  two  columns.  The 
frequency  is  wave  frequency  and  the  PSD  is  normalized  by  the 
resolution  bandwidth  in  radians.  The  last  two  columns  are 
the  data  transformed  to  encounter  frequency  spectra  and 
normalized  by  the  resolution  bandwidth  in  Hertz.  The  encounter 
direction  was  25  off  the  starboard  bow. 
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TABLE  4 


NOAA  WAVE  BUOY 
11-10-78 
COPPER  HARBOR 


Hertz  PS^  Radian  dfo  Hertz  PSIJ  Hertz 

Fo  cni  -sec  /  e  pe  cm  _sec 


0.016 

17.841 

.883 

.017 

98.98 

0.033 

18.581 

.796 

.037 

92.93 

0.049 

22.573 

.  724 

.058 

102.68 

0.065 

25.115 

.  664 

.081 

104.78 

0.081 

22.256 

.613 

.106 

85.72 

0.098 

20.421 

.  567 

.135 

72.75 

0.114 

21.219 

.  530 

.165 

70.66 

0.130 

23.327 

.497 

.196 

78.84 

0.146 

27.276 

.  468 

.229 

80.21 

0.163 

36.897 

.441 

.266 

102.24 

0.179 

53.713 

.418 

.304 

141.07 

0.195 

71.895 

.397 

.343 

179.34 

0.211 

78.833 

.  378 

.  384 

187.23 

0.288 

72. 141 

.360 

.430 

163.18 

0.244 

67.508 

.  345 

.476 

146.34 

0.260 

|  74.567 

.  331 

.523 

155.08 

0.276 

87.490 

.  318 

.573 

174.79 

0.293 

99.687 

.305 

.627 

191.04 

0.309 

99.477 

.294 

.681 

183.76 

0.325 

92.602 

.283 

.736 

164.66 

0.341 

67.492 

.274 

.794 

116.19 

0.358 

62.036 

.264 

.857 

102.90 

0.374 

57.648 

.256 

.919 

92.73 

0.390 

54.138 

.248 

.982 

84.36 

0.406 

53.048 

.240 

1.048 

80.00 

0.423 

51.084 

.233 

1.120 

74.79 

0.439 

48.446 

.226 

1.189 

68.79 

0.455 

48.698 

.220 

1.261 

67.32 

0.471 

52.240 

.214 

1.335 

70.24 

0.489 

54.485 

.208 

1.420 

71.21 

Power  Spectral  Densities  expressed  in  terms  of  encounter 
frequencies  (fe)  for  runs  6,  7,  10,  14,  23,  27,  29,  and  30 
are  shown  in  figures  11  through  18.  Encounter  directions 
(c.f.  eg.  (8))  are  taken  to  be  the  difference  between  the 
ship's  heading  and  the  observed  wave  direction.  Note  that 
the  ship’s  resonant  springing  frequency  is  0.32  Hz  and  the 
critical  encounter  frequencies  are  of  the  order  of  1.0  or 
1.5  Hz  for  the  Collins  or  NRL  radars. 


CONCLUSIONS 

The  evaluation  of  the  wave  measurement  system  was  done 
by  a  comparison  of  the  significant  wave  heights  derived  from 
the  radar  measurements  and  visual  estimates  as  listed  in 
Table  3.  In  addition,  the  spectra  of  the  radar  range 
measurements  and  the  vertical  accelerometer  show  a  resonance 
of  the  vessel  at  the  same  0.32  Hz  frequency  as  an  earlier 
study  (reference  1) .  A  good  estimate  of  the  wave  encounter 
frequency  spectra  and  SWH  measurements  are  obtained  by  the 
wave  measurement  system  as  indicated  by  the  results. 


The  wave  measurement  system,  as  described  in  this 
report,  aboard  the  M/V  STEWART  J.  CORT  will  be  used  by  the 
Naval  Ship  Research  and  Development  Center  (NSRDC)  for  a 
study  of  the  response  of  the  vessel  to  waves  encountered. 

The  field  measurements  by  NSRDC  will  consist  of  simultaneous 
measurements  of  hull  pressure,  structural  stress  and  wave 
data. 
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APPENDIX  B 


TABLE  B-l 
SERIAL  NO.  6 

GREAT  LAKES  WAVE  HEIGHT  PROJECT 

DATE:  11  Oct  78  TIME:  19:06:05 


DURATION  OF  RUN  IN  MINUTES  IS 
NORTH  LATITUDE  (DD  MM) 

WEST  LONGITUDE  (DD  MM) 
VESSEL'S  SPEED  (MPH  -  XX. X) 
VESSEL'S  HEADING  (DEGREES) 
VESSEL'S  DRAFT  (FEET) 

WIND  DIRECTION  (DEGREES) 

WIND  SPEED  (KNOTS) 

WAVE  DIRECTION  (DEGREES) 

WAVE  HEIGHT  (FEET) 

REMARK  LINE  -  UP  TO  60  CHARACTERS 


20 

47  40 
87  50 
15.3 
100 
0 

180 

16 

195 

5 


DRAFT  FORWARD2 7- 2 ,  MID  27-3,  AFT  27-4 


REMARK  LINE  -  UP  TO  60  CHARACTERS 
SL-7  RADAR  NOT  WORKING,  10.5  MILES  OFF  COPPER  HARBOR 


REMARK  LINE  -  UP  TO  60  CHARACTERS 
PROG  INT  2  TIMES  7:25  &  7:30  DARK,  WAVES  OBS  GUESSES 


POWER  SPECTRAL  DENSITY  SN  #6  11  OCTOBER  1978 


33S-WD 

'(b)  t 


- 


FREQUENCY  OF  ENCOUNTER  -  HERTZ 

FIGURE  11 


TABLE  B-2 


SERIAL  NO.  7 


GREAT  LAKES  WAVE  HEIGHT  PROJECT 
DATE:  11  Oct  78  TIME:  19:39:31 

DURATION  OF  RUN  IN  MINUTES  IS 
NORTH  LATITUDE  (DD  MM) 

WEST  LONGITUDE  (DD  MM) 

VESSEL'S  SPEED  (MPH  -  XX. X) 

VESSEL'S  HEADING  (DEGREES) 

VESSEL'S  DRAFT  (FEET) 

WIND  DIRECTION  (DEGREES) 

WIND  SPEED  (KNOTS) 

WAVE  DIRECTION  (DEGREES) 

WAVE  HEIGHT  (FEET) 

REMARK  LINE  -  UP  TO  60  CHARACTERS 
DRAFT  LOADED 


15 

47  50 
87  50 
15.3 
100 
0 

180 

16 
195 

5 


REMARK  LINE  -  UP  TO  60  CHARACTERS 
SL-7  RADAR  NOT  WORKING 

REMARK  LINE  -  UP  TO  60  CHARACTERS 
DARK  WAVE  OBS  JUST  GUESSES  PROB  6  7  8  FT  WAVES 


3Q(«) 


. . 


'  ■  : *\  • 


— 


POWER  SPECTRAL  DENSITY  SN  #  7  11  OCTOBER  1978 

_ 1939  HOURS 


FREQUENCY  OF  ENCOUNTER  -  HERTZ 

FIGURE  12 


TABLE  B-3 


SERIAL  NO.  10 


GREAT  LAKES  WAVE  HEIGHT  PROJECT 
DATE:  12  OCT  78  TIME:  15:24:19 

DURATION  OF  RUN  IN  MINUTES  IS 
NORTH  LATITUDE  (DD  MM) 

WEST  LONGITUDE  (DD  MM) 

VESSEL'S  SPEED  (MPH  -  XX. X) 

VESSEL'S  HEADING  (DEGREES) 

VESSEL'S  DRAFT  (FEET) 

WIND  DIRECTION  (DEGREES) 

WIND  SPEED  (KNOTS) 

WAVE  DIRECTION  (DEGREES) 

WAVE  HEIGHT  (FEET) 


20 

45  52 
84  9 

14.5 
264 
0 

270 

40 

260 

6 


REMARK  LINE  -  UP  TO  60  CHARACTERS 
1.6  MILES  OFF  MARTIN  REEF 


REMARK  LINE  -  UP  TO  60  CHARACTERS 
DRAFT,  LOADED,  NO  SL-7  DATA 

REMARK  LINE  -  UP  TO  60  CHARACTERS 
HEAD  SEAS,  SPRAY  OVER  BOW 


POWER  SPECTRAL  DENSITY  SN  #  10  12  OCTOBER  1978 

1529  HOURS 


FREQUENCY  OF  ENCOUNTER  -  HERTZ 

FIGURE  13 


TABLE  B-4 


SERIAL  NO.  14 


GREAT  LAKES  WAVE  HEIGHT  PROJECT 
DATE:  22  OCT  78  TIME: 

DURATION  OF  RUN  IN  MINUTES  IS 
NORTH  LATITUDE  (DD  MM) 

WEST  LONGITUDE  (DD  MM) 

VESSEL'S  SPEED  (MPH  -  XX. X) 

VESSEL'S  HEADING  (DEGREES) 

VESSEL'S  DRAFT  (FEET) 

WIND  DIRECTION  (DEGREES) 

WIND  SPEED  (KNOTS) 

WAVE  DIRECTION  (DEGREES) 

WAVE  HEIGHT  (FEET) 


14:50:33 

20 

46  50 
85  12 
15.4 
292 
17 
320 
35 
310 
8 


REMARK  LINE  -  UP  TO  60  CHARACTERS 
SL-7  IN  SMOOTH  MODE,  BALLAST  CODE  3,  7  MILES  CRISP  PT 


REMARK  LINE  -  UP  TO  60  CHARACTERS 
DRAFT  FWD  17-6,  MID  19,  AFT  21 


REMARK  LINE  -  UP  TO  60  CHARACTERS 
SOME  VESSEL  MOTION,  SPRAY  OVER  SIDE  AT  TIMESDEEPTH  233  FT 
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POWER  SPECTRAL  DENSITY  SN  #  W  22  OCTOBER  1978 


FREQUENCY  OF  ENCOUNTER  -  HERTZ 

FIGURE  14 


TABLE  B-5 


SERIAL  NO.  23 

21:24:22 

5 

47  32 
87  40 
1.6 
287 
19 
65 
16 
330 
4 


REMARK  LINE  -  UP  TO  60  CHARACTERS 
7  MILES  EAST  OP  COPPER  HARBOR 

REMARK  LINE  -  UP  TO  60  CHARACTERS 
10  MIN  TILL  COURSE  CHANGE 


GREAT  LAKES  WAVE  HEIGHT  PROJECT 
DATE:  09  NOV  78  TIME: 

DURATION  OF  RUN  IN  MINUTES  IS 
NORTH  LATITUDE  (DD  MM) 

WEST  LONGITUDE  (DD  MM) 

VESSEL'S  SPEED  (MPH  -  XX. X) 

VESSEL'S  HEADING  (DEGREES) 

VESSEL'S  DRAFT  (FEET) 

WIND  DIRECTION  (DEGREES) 

WIND  SPEED  (KNOTS) 

WAVE  DIRECTION  (DEGREES) 

WAVE  HEIGHT  (FEET) 

REMARK  LINE  -  UP  TO  60  CHARACTERS 
CODE  4  BALLAST 
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POWER  SPECTRAL  DENSITY  SN  U  23  9  NOVEMBER  1978 

2102  HOURS 


FREQUENCY  OF  ENCOUNTER  -  HERTZ 

FIGURE  15 


TABLE  B-6 


SPECTRAL  DENSITY  SM  ft  27  12  NOVEMBER  1978 

11:25  AND  11:33 


FREQUENCY  OF  ENCOUNTER  -  HERTZ 

FIGURE  16 


TABLE  B-7 


SERIAL  NO.  29 


GREAT  LAKES  WAVE  HEIGHT  PROJECT 

DATE:  16  NOV  78  TIME:  22:52:33 


DURATION  OF  RUN  IN  MINUTES  IS 
NORTH  LATITUDE  (DD  MM) 

WEST  LONGITUDE  (DD  MM) 
VESSEL'S  SPEED  (MPH  -  XX. X) 
VESSEL'S  HEADING  (DEGREES) 
VESSEL'S  DRAFT  (FEET) 

WIND  DIRECTION  (DEGREES) 

WIND  SPEED  (KNOTS) 

WAVE  DIRECTION  (DEGREES) 

WAVE  HEIGHT  (FEET) 

REMARK  LINE  -  UP  TO  60  CHARACTERS 


4 

47  29 
87  54 
16.3 
270 
17 
315 
20 
315 
4 


COPPER  HARBOR  1.5  MILES,  CODE  3  F  17,  M  18-11,  AFT  21-3 


REMARK  LINE  -  UP  TO  60  CHARACTERS 
JUST  WEST  OF  BUOY 


REMARK  LINE  -  UP  TO  60  CHARACTERS 
CHANGED  HEAD  TO  265 


75  (*  ) 


POWER  SPECTRAL  DENSITY  SN  ft  29  16  NOVEMBER  1978 

221?  HOURS 


figure  17  FREQUENCY  OF  ENCOUNTER  -  HE&TZ 


TABLE  B-8 


SERIAL  NO.  30 


GREAT  LAKES  WAVE  HEIGHT  PROJECT 
DATE:  23  NOV  78  TIME 

DURATION  OF  RUN  IN  MINUTES  IS 
NORTH  LATITUDE  (DD  MM) 

WEST  LONGITUDE  (DD  MM) 

VESSEL’S  SPEED  (MPH  -  XX. X) 
VESSEL'S  HEADING  (DEGREES) 

VESSEL'S  DRAFT  (FEET) 

WIND  DIRECTION  (DEGREES) 

WIND  SPEED  (KNOTS) 

WAVE  DIRECTION  (DEGREES) 

WAVE  HEIGHT  (FEET) 


23:46:49 

15 

41  44 
87  11 
15.3 
343 
17 
310 
40 
290 
6 


REMARK  LINE  -  UP  TO  60  CHARACTERS 
BALLAST  CODE  3,  DEP  72  FT,  STRIP  CHART  OF  BOTH  RADARS 


REMARK  LINE  -  UP  TO  60  CHARACTERS 
LEFT  BURNS  HAR,  GOING  UP  WEST  SIDE  TO  STAY  IN  LEE 

REMARK  LINE  -  UP  TO  60  CHARACTERS 
SPRAY,  SHIP  MOTIONS,  WAVE  IMPACT,  15  GOOD  MIN  OF  DATA 
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APPENDIX  A 
TABLE  A-l 


SUMMARY  OF  OPERATIONS  ABOARD  THE  VESSEL 
March  1978 

*  Installation  of  System 


April  1978 

*  Installation  of  cable  connectors 

*  Initial  check  of  system  on  vessel 


June  1978 

*  Return  computer  to  Washington,  D.C.  PHASE  I 


July  &  August  1978 

*  Calibrate  System  while  in  drydock 

*  Attach  longer  boom  for  bow  radar 


October  &  November  1978 

*  Collect  data 

April  1979 

*  Install  longer  boom  for  Collins  Radar 

May  &  June  1979  PHASE  II 

*  Trial  run 

October  1979 

*  Repair  NRL  Radar  Transmitter 
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APPENDIX  B 


LOG  SHEETS  FOR  REDUCED  DATA 


(Reduced  Data  Sheets  were  transferred 
to  their  associated  PSD  plot.  Page 
numbers  39  -  46  is  maintained  here 
for  continuity) 

HEW 

P_age  # 

•29(a) 

.30(a) 

•31(a) 

•32(a) 

•33(a) 

•34(a) 

•35(a) 

•36(a) 


OLD 

page  # 

39.  .  . 

40.  .  . 

41.  .  . 

42.  .  . 

43.  .  . 

44.  .  . 

45.  .  .  . 
46.... 


APPENDIX  C 


NRL  RADAR  SYSTEM 


The  NRL  radar's  high  resolution  of  30  centimeters  is 
obtained  by  using  a  very  narrow  pulse  of  2  nanoseconds. 
Figure  C-l  is  a  block  diagram  of  the  radar  system.  The 
important  system  characteristics  are  listed  in  Table  C-l. 

TABLE  C-l 


Wavelength 
Pulse  Width 

Peak  Transmitted  Power 
Pulse  Repetition  Rate 
Antenna  Diameter 
Receiver  Noise  Figure 
Equivalent  Pulse  Processing  Rate 


3  Centimeters 
2  Nanoseconds 
100  Watts 
8000  per  second 
60  Centimeters 
7  db 

80  per  second 


The  high  frequency  components  for  the  transmitter  and 
receiver  are  located  on  the  antenna  pedestal  on  the  ship's 
deck  about  21  meters  above  the  water  line.  The  antenna  is 
pointed  ahead  and  tilted  down  and  out  about  25  degrees  with 
respect  to  the  vertical. 


The  control  and  display  circuits  are  located  inside  the 
pilot  house  in  a  standard  half  rack.  The  timing  and  control 
signals  are  derived  by  using  both  digital  and  analog  computer 
circuits. 

A  sampling  scope  transforms  the  received  periodic  signals 
which  have  a  bandwidth  of  1  GHz  to  equivalent  time  signals  with 
a  bandwidth  of  8  KHz. 


SAMPLING  SCOPE 

A  detailed  description  of  the  operation  of  the  sampling 
scope  is  found  in  references  6,  7,  and  8.  A  sampling  scope 
requires  a  number  of  repetitions  of  a  periodic  signal  in  order 
to  construct  an  "A"  scope  sweep  display  of  a  single  echo.  In 
the  NRL  system  the  sampling  scope  sweep  range  window  represents 
100  nanoseconds  which  is  equivalent  to  15  meters.  The  signal 
amplitude  from  the  first  pulse  returned  is  sampled  at  the  range 
of  the  first  nanosecond  of  the  sweep  window.  The  second  pulse 
is  sampled  at  the  second  nanosecond  and  so  forth  for  100  return 
pulses,  which  together  constitute  the  scope  trace.  The  scope 
sweep,  therefore,  represents  a  composite  picture  of  a  single 
echo.  The  transmitter  repetition  rate  of  8000  pulses  per 
second  allows  the  reproduction  of  80  high  resolution  echoes 
per  second. 
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SAMPLING 

SCOPE 


RADAR  SYSTEM 


AUTOMATIC  GAIN  CONTROL 


The  block  diagram  for  the  automatic  gain  control  is  shown 
in  Figure  C-2.  The  amplitude  of  the  return  has  rapid  changes 
of  about  100  Hz.  The  average  amplitude  also  changes  because  of 
the  changing  aspect  angle  of  the  antenna  beam  and  the  water 
surface.  This  changes  the  amplitude  of  average  return  signal 
by  about  ten  to  one.  The  AGC  circuit  controls  the  RF  gain  of 
the  receiver  through  a  voltage  controlled  pin  diode  attenuator 
and  can  compensate  for  average  amplitude  changes  of  about  10 
Hz,  but  the  more  rapid  fading  of  the  pulse  is  not  corrected  by 
the  AGC. 

RANGE  TRACKER 

The  block  diagram  for  the  range  tracker  is  shown  in  Figure 
C-3.  The  range  tracker  uses  the  programmable  feature  of  the 
Tektronix  3T5  sampling  sweep  unit  (reference  8) .  The  digital 
control  for  the  sweep  range  window  is  a  counter  which  counts 
from  zero  for  the  beginning  of  the  sweep  up  to  one-hundred  at 
the  end  of  the  sweep.  The  range  position  of  the  start  of  the 
sweep  range  window  is  controlled  by  a  programmable  delay  circuit 
in  the  scope.  This  delay  is  either  increased  or  decreased 
depending  on  the  position  of  the  leading  edge  of  the  return 
pulse  within  the  sweep.  The  position  of  the  leading  edge  is 
detected  by  a  peak  amplitude  discriminator  circuit.  The 
discriminator  produces  a  range  pulse  at  the  time  of  the  echo 
pulse  leading  edge.  This  range  pulse  stores  the  digital  number 
in  the  range  counter.  This  number  is  then  used  to  increment  a 
bi-directional  range  delay  counter  which  in  turn  is  used  to 
program  the  sweep  delay  of  the  sampling  scope.  This  variable 
delay  keeps  the  echo  centered  in  the  range  window. 

The  decision  to  increase  or  decrease  the  range  delay  is 
accomplished  by  utilizing  a  digital  range  comparator.  The 
comparator  is  strobed  at  the  end  of  every  sweep.  The  delay 
is  decreased  by  one  nanosecond  if  the  range  pulse  comes  before 
the  range  counter  reaches  55.  The  delay  is  increased  one 
nanosecond  if  the  range  pulse  comes  after  the  counter  reaches 
55.  If  no  range  pulse  is  detected  during  a  sweep,  the  delay 
is  increased  by  two  nanoseconds  which  allows  the  radar  system 
to  search  through  an  extended  range  for  the  echo. 

CONCLUSION 

This  radar  system  measures  the  range  from  the  ship  to  the 
ocean  surface.  The  range  is  digitized  in  increments  of  one 
nanosecond  which  is  approximately  15  centimeters,  however,  pulse 
jitter  and  the  uncertainty  of  locating  the  leading  edge  of  the 
echo  reduces  the  range  accuracy  to  about  30  centimeters.  Infor¬ 
mation  about  the  ocean  wave  heights  can  be  obtained  by  subtract¬ 
ing  the  independently  measured  range  changes  due  to  the  ship's 
motion  from  the  radar  data. 
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A.G.C.  SIGNAL  TO  RECEIVER 


AUTOMATIC  GAIN  CONTROL 


TRIGGER  TO  TRANSMITTER  ECHO  RETURN 


NRL  RADAR  SYSTEM  MODIFICATIONS 

The  NRL  microwave  radar,  originally  designed  to  measure 
ocean  waves  from  a  SL-7  class  containership,  was  modified 
for  operations  aboard  the  1000  foot  Great  Lakes  ore  vessel. 

The  modifications  provide  for  the  measurement  of  the  total 
distance  to  the  surface  and  also  to  convert  the  data  to  a 
suitable  rate  in  digital  format  for  recording  by  the  PDP-11/03 
computer. 

The  original  radar  system  used  a  10,000  per  second  timing 
circuit.  This  was  replaced  by  a  more  precise  timing  circuit 
that  uses  a  crystal  oscillator.  This  provides  transmitter 
pulses  at  the  rate  of  10,560  per  second.  One  hundred  pulses 
are  used  for  each  trace  of  the  sampling  scope  and  thirty-two 
are  used  as  time  for  the  sampling  scope  to  reset.  Thus, 
there  are  eighty  pulse  returns  received  and  processed  per 
second.  Eight  successive  range  estimates  are  averaged  to 
produce  ten  samples  of  radar  range  per  second,  each  with  a 
precision  of  1.875  centimeters. 

The  original  SL-7  version  of  the  NRL  microwave  radar 
system  used  a  binary  counter  followed  by  a  circuit  to  convert 
the  digital  value  to  an  analog  voltage  proportional  to 
changes  in  distance  to  the  surface.  The  data  recording  of 
the  Great  Lakes  wave  measurement  system  uses  a  digital 
format  to  match  the  digital  computer  used.  This  was  ac¬ 
complished  by  replacing  the  relative  range  counter  with  two 
digital  circuits.  The  first  circuit  converted  the  binary 
coded  decimal  form  of  the  distance  data  to  straight  binary. 

The  second  digital  circuit  summed  eight  successive  range 
values  and  transmitted  them  to  the  onboard  PDP-11/03  com¬ 
puter  via  its  digital  input  board  type  DRV  11. 

The  revised  assignments  for  the  location  of  the  circuit 
boards  in  the  Great  Lakes  version  of  the  NRL  microwave  radar 
system  is  shown  in  Figure  C-4.  Slot  7  contains  the  new 
crystal  oscillator  timing  circuit.  Slot  21  now  contains  the 
digital  circuit  board  that  converts  the  binary  coded  decimal 
range  values  to  straight  binary.  Slot  22  holds  the  digital 
circuit  board  that  sums  eight  successive  range  values  and 
transmits  them  to  the  PDP-11/03  computer. 

The  wiring  diagram  for  the  new  timing  circuit  card  is 
shown  in  Figure  C-5.  The  LS  160 's  are  decade  dividers.  The 
LS  30  multiple  input  gate  provides  an  output  pulse  for 
every  947  input  pulses.  The  LS  265  provide  four  parallel 
outputs  to  the  other  circuits  in  the  radar  system. 

The  binary  counter  circuit  boards  and  the  associated 
wiring  that  were  located  in  slots  21,  22,  and  27  were  removed. 
The  circuit  board  that  converts  the  binary  coded  decimal 
(BCD)  range  data  to  straight  binary  form  is  now  in  slot  21. 
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CARD  LAYOUT,  NRL  RADAR 
Fig.  C-4 


The  wiring  diagram  for  this  circuit  board  is  shown  in  Figure 
C-6 .  The  twelve  lines  that  represent  the  BCD  numbers  are 
inputs  on  the  left  of  Figure  C-6.  The  straight  binary 
representation  for  the  range  appears  as  the  ten  wires  on  the 
right  side  of  Figure  C-6.  The  numbers  inside  the  little 
squares  refer  to  pin  numbers  on  the  circuit  card.  The  least 
significant  digit  is  at  the  top  of  Figure  C-6  for  both  the 
BCD  input  and  the  binary  output.  The  counter  and  logic 
gates  at  the  bottom  of  Figure  C-6  provide  for  an  exact  time 
between  sampling  scope  data  sweeps  of  32  timing  pulses. 

This  time  between  sweeps  is  used  to  reset  the  sampling  scope 
and  to  update  the  range  counter. 

The  wiring  diagram  for  the  eight-word  summer  is  shown 
in  Figure  C-7.  The  three  LS  283' s  are  four  binary  digit 
summers  with  carry. 

The  intermediate  results  are  stored  in  circuits  LS  175 
just  below  the  summers.  Thus,  the  intermediate  values  are 
added  to  the  current  range  value.  After  the  eighth  sum  the 
value  is  transferred  to  the  output  storage  registers  at  the 
bottom  of  the  figure.  These  twelve  data  outputs  along  with 
the  data  ready  pulse  are  transmitted  directly  to  the  PDP- 
11/03  computer  via  the  system  cable  VI.  The  cable  connector 
is  a  25  pin  cannon  type-D  located  on  the  top  of  the  radar 
system  electronic  chassis. 
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APPENDIX  D 


WAVE  MEASUREMENT  SYSTEM  INTERCONNECTIONS 

The  electrical  cables  that  connect  the  principal  com¬ 
ponents  of  the  wave  measurement  system  with  the  PDP-11/03 
computer  are  shown  in  Figure  D-l.  The  NRL  microwave  radar 
chassis,  the  signal  distribution  chassis,  and  the  PDP-11/03 
computer  are  located  in  the  pilot  house.  The  NRL  microwave 
radar  transceiver  and  a  pair  of  accelerometers  are  mounted 
on  top  of  the  pilot  house.  The  Collins  radar  horns  and 
another  pair  of  accelerometers  are  mounted  on  a  4.5  meter 
long  boom  over  the  bow  at  the  main  deck.  The  Collins  radar 
and  electronic  circuits  are  contained  in  a  weather  proof 
chassis  at  the  bow.  Tables  D-l  through  D-7  identify  cable 
signal  allocations. 

Figure  D-2  is  the  wiring  -’ithin  the  signal  distribution 
chassis  in  the  pilot  house.  The  digital  values  of  NRL  radar 
range  are  on  jack  J-l.  These  are  directed  to  jack  J-2  and 
the  computer  when  the  test-operate  switch  is  in  the  operate 
position.  The  analog  voltage  signals  as  defined  in  Table  D-7 
are  directed  to  the  computer  via  the  multiple  pole  T-bar 
switch  when  the  test-operate  switch  is  in  the  operate  position. 

The  test  generator  contains  a  digital  counter  that 
increases  one  count  ten  times  a  second.  The  digital  form  of 
the  counter  is  recorded  by  the  computer  instead  of  the  NRL 
radar  signal  when  the  test-operate  switch  is  in  the  test 
position.  The  digital  value  is  also  converted  to  an  analog 
voltage  and  directed  to  all  the  analog  inputs  of  the  computer 
when  the  test-operate  switch  is  in  the  test  position. 

Normally  the  test-operate  switch  is  always  in  the  operate 
position.  If  a  test  of  the  computer  recording  is  needed, 
the  switch  is  placed  in  the  test  position  and  the  computer 
should  record  via  the  DRV-11  input  a  count  that  increased  by 
one  each  sample  until  the  full  scale  value  is  reached,  at 
which  time  the  count  resets  the  zero  and  starts  over  again. 

The  analog  inputs  on  the  ADV-11  should  follow  in  a  similar 
stair  step  fashion. 

NOTE:  PDP-11/03  Real-Time  I/O  OPTIONS 

DRV-11  :  Parallel  line  interface  unit.  16-bit  diode  input;  16-bit 
latched  drive  output.  Protocol  &  control  signal.  Board  size 
21.6  cm  X  12.7  cm.  Mounting  Code:  Single  ll/o3  slot. 

ADV-11  :  12-bit,  16  channel  single  ended  or  8  channel 
quasi-differential  analog  -  digital  converter. 
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RADAR 


TKST  I'AMil  wm% 


SFF  TABLF  D-8 


CABF T  f  ROM  DRV  II  OUTI’IIT 


TABLE  D-l 


CABLE  I 

A  +6VDC 

B  Ground 

C  Bow  vertical  accelerometer  -  Low 

D  Bow  vertical  accelerometer  -  High 

E  Bow  horizontal  accelerometer  -  Low 

F  Bow  horizontal  accelerometer  -  High 


TABLE  D-2 

CABLE  II  a,  b 

These  are  two  co-ax  cables  that  connect  the  Collins  radar 
to  the  transmitter  antenna  and  to  the  receiver  antenna- 


TABLE  D-3 


CABLE  III 

A  &  R  115  vac 

B  &  P  115  vac 

F  Collins  low  return  indicator 

G  Chassis  GND. 

J  Collins  radar  range 

K  Chassis  GND. 

U  Bow  vertical  accelerometer  -  High 

V  Bow  vertical  accelerometer  -  Low 

c  Bow  horizontal  accelerometer  -  High 

b  Bow  horizontal  accelerometer  -  Low 


61 


TABLE  D-4 


Cable  IV 

PIN 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


22  (Co-Ax) 

23  (Co-Ax) 

24  (Co-Ax) 


SIGNAL 

Spare  Ground 
Shield  Ground 

Transmitter  Power,  +700  VDC 

Shield  Ground 

Therostat 

Transmitter  Power,  +280  VDC 
+12  VDC 

Transmitter  Power,  -5  VDC 
Receiver  Gain  Control  Return,  Ground 
Starboard,  SPARE 
Port,  SPARE 

Low  Voltage  Power  Ground 

+12  VDC  Sense 

Phone  Ground,  SPARE 

Transmitter  Power  Ground 

Receiver  Gain  Control 

Phone,  SPARE 

110  VAC  Common 

-12  VDC  Sense 

-12  VDC 

110  VAC  Hot 

Transmitter  Trigger  Pulse 
Delayed  Scope  Trigger 
Video  Signal  Return 


TABLE  D-5 


CABLE  V 

PIN 

A 

B 

D 

R 

F 

G 


SIGNAL 

Chassis  Ground 
Horizontal  Accelerometer 

Vertical  Accelerometer 
Chassis  Ground 

110  VAC 
110  VAC 
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TABLE  D-6 


CABLE  VI 
25  PIN 

CANNON  TYPE  'D' 


PIN 

SIGNAL 

1 

Bo 

2 

B1 

3 

B2 

4 

B3 

5 

B4 

6 

B5 

7 

B6 

8 

B7 

9 

B8 

10 

B9 

11 

B10 

12 

B11 

13 

"Three  Missing  Range  Pulses 

14 

"50  Missing  Range  Pulses" 

15 

SPARE 

16 

SPARE 

17 

SPARE 

18 

SPARE 

19 

DATA  READY 

20 

DATA  TRANSMITTED 

21 

SPARE 

22 

SPARE 

23 

SPARE 

24 

SPARE 

25 

GROUND 
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TABLE  D-7 


J-3  and  J-4 


PIN  SIGNAL 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 
17 


Bow  vertical  accelerometer  -  High 
Bow  vertical  accelerometer  -  Low 
Bow  vertical  accelerometer  -  High 
Bow  horizontal  accelerometer  -  Low 
SPARE 
SPARE 

Pilot  House  Top  vertical  accelerometer  -  High 
Pilot  House  Top  vertical  accelerometer  -  Low 
Pilot  House  Top  horizontal  accelerometer  -  High 
Pilot  House  Top  horizontal  accelerometer  -  Low 
Collins  Radar  Range  Signal  -  High 
Collins  Radar  Range  Signal  -  Low 
Not  Used 

Collins  Radar  Signal  Strength  -  High 
Collins  Radar  Signal  Strength  -  Low 
SPARE 
SPARE 


ACCELEROMETER  ELECTRONICS 

The  accelerometer  pair  on  top  of  the  pilot  house  is  in 
a  weather  proof  chassis  attached  to  the  NRL  microwave  radar 
transceiver  stand.  The  electronics  used  to  condition  the 
accelerometer  outputs  are  also  inside  the  same  chassis.  The 
wiring  diagram  for  the  amplifiers  are  shown  in  Figure  D-3. 

The  vertical  accelerometer  signal  is  balanced  by  the  adjustment 
of  the  poteniometer  on  U5  of  the  board.  Thus,  the  vertical 
accelerometer  voltage  recorded  is  near  zero  volts  for  the 
stationary  condition  of  the  vessel. 

The  electronic  circuits  that  condition  the  bow  mounted 
accelerometers  are  located  in  the  weather  proof  chassis  at 
the  bow  that  also  contains  the  Collins  Radar.  Figure  D-4  is 
the  wiring  diagram  of  the  conditioning  circuits  for  the  bow 
mounted  accelerometer  signals.  The  accelerometers  are  located 
in  a  weather  proof  container  at  the  end  of  the  4.5  meter  bow 
boom.  The  vertical  accelerometer  signal  is  balanced  by  adjust¬ 
ment  of  the  poteniometer  on  44  for  a  value  near  zero  volts  for 
the  stationary  condition  of  the  vessel. 
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BOW  BOOM  ACCELEROMETER  SIGNAL  CONDITIONING  WIRING 


APPENDIX  E 


COLLINS  RADAR 

The  Collins  Radar  was  designed  as  an  aircraft  radio 
altimeter.  It  is  a  frequency  modulated  transmitted-receiver 
that  operates  in  the  4250  to  4350  megahertz  frequency  range. 
The  radar  beams  the  signal  to  the  water  surface,  receives 
the  reflected  signal  and  converts  it  to  an  analog  voltage 
proportional  to  the  distance  to  the  surface. 

The  Collins  Radar  uses  separate  antennas  for  transmis¬ 
sions  and  reception  of  signals.  The  standard  antennas  for 
the  Collins  aircraft  installation  have  been  replaced  by 
larger  horn  antennas.  These  provide  a  narrow  beam  width  and 
therefore  achieve  better  spatial  resolution  as  discussed 
earlier  in  this  report. 


TABLE  E-l 


Collins  System  Specifications 


Transmitter  output 
Power 

Center  Frequency 
Modulation 

FM  Deviation  Peak-to-peak 

Antenna  Gain 
Antenna  Beam  Width 


150  milliwatts 
4300+  15  MHz 


98.4  MHz 
20  db 

20  degrees 


Svstron  Donner  Antenna  Horns 
Model  HBK- 520-20-1 
Center  frequency  4.3  GHz 
BW  100  MHz 

Gain  -  nominal  20  db 


COLLINS  RADIO  CO. 

Radar  Altimeter  (860F-2) 
ALT- 50 
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Unit  Weight  5.2  pounds  (2.34  ky) 


COLLINS  RADAR  WIRING 


The  wiring  for  the  Collins  Radar  power  and  signal 
conditioning  circuits  are  shown  in  Figure  E-l.  The  offset 
potentiometer  at  the  top  of  the  figure  is  adjusted  to  produce 
an  output  to  the  computer  near  zero  volts  for  the  antenna 
■horns  mounted  at  the  end  of  the  bow  shown  viewinq  the  sur¬ 
face  while  in  port. 
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COLLINS  RADAR  SIGNAL  CONDITIONING  WIRING 


APPENDIX  F 


FLOW  DIAGRAM  for  M/V  CORT  Data  Analysli  ProgratR 


START 


INPl'T  PARAMETERS 
OK  ANALYSIS  .) 


0 


YES 


•'  PRINT  DATA 


tOMPlTE  AND  PRIM 
VARIANCE 


1 


NO  OF  PTS 

APPLY 

PT  SPACING 

- 

\  B-H  WINDOW? 

i 

x 

WINDOW 
!  DATA 


I 


COMPUTE 
!  FFT 


T“J, _ 

'  PRINT 
i  COEFF 

i 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


CORT  0 A T ft  ANALYSIS  PROGRAM  FOR  THE  79-79  MEASUREMENT  SEASON 
WRITTEN  BY  DAVID  A.  WALOEN 

OFFICE  OF  RESEARCH  AND  DEVELORMENT 
U.  S.  COAST  GUARC 


note:  all  calculations  and  results  are  IN  COMPJTER  COU\ItS 

CI.E.  THERE  IS  NO  CONVERSION  TO  ENGINEERING  UNITS) 


c 

DIMENSION  INDIG<  256) • INAL G! 8.256 ),X!2099>  #Y I2C9 7) 

DIMENSION  INPUT!  2304) »  UL I  22 )  »NLIM!22»9) 

DIMENSION  1 0 ATE!  9 ) » I T IHE ! 9  >  » I L A T ! 2 ) » ILONG! 2 ) t  L  K'JP ! 10) 

E DU I  VALENCE  ! HEADER, IDATE)  .  ( HEADER !6  > ♦ITIME)»1INAME«INAV> 
EQUIVALENCE  !  HEADER  !  1 0  )  »  I  OUR  A  T)  ,  1  H£  AOER  I  1 1  ),»  I L  A  T  ) 

EQUIVALENCE  !  HE  AO  rR ! 1 3 )  ,  I LCN 3 ) , 1  HE ADER ! 1 5 > , S  =  E E 0 > 

EQUIVALENCE  ! HEADER  1 1 7 ) , I  HE A D ) , ! HE  AO ER  1  1  9 ) *  ID R A ^T ) 

EQUIVALENCE  !  HEAD™!  19),  IUNOIR),  !  HEADER!  20)  »IWN03P) 

'  EQUIVALENCE  !  HEAD  ER  !  2  i  )  ,  I U  AV  OR  )  t  <  HE  A  OCR  !  22  )  ,  I  .'A  VH  T ) 

EQUIVALENCE  ! INPUT, I  NO  I G) t ( I N»U T !  2 57 ) , I N AL G > 

L  03 1  C  AL  •  1  NA  M  A  !  9  ) 

COHMON/BLKl/INPUT 
BYTE  IDATE.ITIME 

INTEGER*?  HEA0ERI256),  INAwr!  40)  *  N  A*»E  1  9  )  ,  BLOCK 

DATA  IN  AT/ JR  0X1 , JR  ROD, JR  AO  ,  J= CAT, 3*0X1, 3RRDD.J°4l  .JRDAT.J’C* 
15RRDO,3RA2  , 3R D A T , 3RDX 1 , 3RR DC. 3 R 4 3  ,  3R D A T , 3 J OX  1 , 3 R R 0 D. 3R A 4  ,3  = 
23R  DX1,3RRD0*3RA5  , JRQ A  T  ,  JRCX1 , JRRDD, 3R  AS  ,  39 D A T  ,  JR DX  l , J R  =  DC » 


J3RA7  , 

3R04  T , JROX1, 3RROC*3RAR  , JRDAT.JRCXI.JRRDO.JRAH 

930 

rQR H  a  T 

(  • 

*  ,  I  6  »5X  ,»I6) 

901 

FORMAT 

(/) 

910 

-CRMAT 

(  t 

•  »10X,*r.PEAT  LAKES  WAVE  HEIGHT  PROJECT'//) 

911 

format 

(  *  ♦ 

‘.l0X,»O4TE*,3X,9Al,9X,»TIME»,3X,3Al,/> 

912 

r  ORM  A  T 

(  « 

• *10X, ‘DURATION  Cc  RUN  IN  MINUTES  IS‘,3*,f> 

913 

c  DR  MAT 

(  • 

'.12) 

919 

FORMAT 

(  9 

•«10X, 'NORTH  LATITUDE  !OD  MM)',9X,J> 

915 

tORMAT 

(  • 

',12, IX, 12) 

916 

-ORMAT 

(  • 

' *  1 0  X  » • Ur  ST  LONGITUDE  IDO  MM)»,9X,$) 

91  7 

-ORMAT 

(  » 

'  » F9 •  1 ) 

919 

rORMAT 

(  t 

• ,10X,2YHVESSEL'S  SPEED  1MRH  -  XX.K),JX,*> 

919 

=  OR  M A  T 

|  t 

' *10X,26HVESSEL*S  HEADING  I  DEGREES )•  5X  ,  t ) 

920 

-ORMA  T 

(  9 

'.I  3) 

921 

“0RM4  T 

(  • 

•,10X,21HVESSEL'S  DRAFT  ! FEET >  ,  1 1 X , S > 

922 

FORMAT 

(  1 

' ,10X,'WI ND  DIRECTION  I  DEGREES ) • , 7X  ,  *  ) 

923 

=  0RMA  T 

(  1 

' • 1  OX, ' WI NO  SPEEO  IKNOTS)'  « 1 9 X , t ) 

929 

FORMAT 

(  1 

• ,10X, 'WAVE  DIRECTION  1  DEGREES )», 7X , * ) 

925 

-ORMAT 

(  • 

' »10X,'WAVE  HEIGHT  IFEET)' ,19X,t> 

926 

FORMAT 

(  t 

REMARK  LINE  -  UP  TO  SO  CHARACTERS') 

92  7 

c0RM A  T 

(  9 

' , 32 A2  ) 

CALL  LOCK 

IFET  : 

IFETCH! INAME) 

IF! IFET 

.NE 

.0)  STOP  'FATAL  ERROR  FETCHING  HANOLER' 

ICHAN=IGETC! ) 

=1=3.19159265 
TYPE  931 

931  eORMAT ! •  WANT  CATALOG  OF  OISK?  TYPE  0  FOR  NO*  1  FOR  YES'*5X,«) 

ACCEPT  932, ISU 
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932  FORMAT C 13  > 

IF(ISW.E9.0)3O  TO  501 
C 

C  .  PRINT  OUT  LISTING  CF  ALL  HEADER  RECOROS  ON  DISK 

C 

CALL  CLOSEC  ( 1  CHAN ) 

00  1  11=1*10 
ICHAN=IGETC( » 

00  21  J=l,4 
JJ  =  < II-l) *4*J 
21  NAMEC J)= INAME( JJ> 

CALL  R50A3C(9,.N4ME,NA'<A> 

TYPE  939, NAHA 

9C9  FORMAT!*  •/////ROC* *• >//*  FILE  NAME=*«941> 

LKUP  (  m=LOOKUP<  ICHAN.NAME) 

TYP'  948,LKUP(II) 

948  rORH A  T<  *  NO  CF  3LCCKS=*,IS> 

I - ( L KUP( 1 1 )  .L T .0 ) GO  r0  501 
TYPE  910 

ICODE=IREAOW<  256,HEA0r9, l , IChAN  ) 

T  Y°F  911.  IOATE.ITIME 
TYPE  912 

TYPE  913,  I  OUR  AT 

TYPE  914 

TYPE  915,  ILAT 

T YPr  916 

TYP"  915,  ILONG 

TYPE  919 

T  Y°E  017,  SPEED 

TYP"  919 

TYP"  J20 .  IHE40 

T  Y°-  921 

TYPE  913,  IORAFT 
TYPE  922 

TYP'-  920,  IWNOIR 
TYP-  923 

TYPE  913,  IWN3SP 
TYPE  924 

TYPE  920,  IWAVDR 
TYPE  925 

TYPE  913,  IWAVHT 
T YPC  926 

TYPE  927,  (HEADER (I). 1=33, 64) 

TYPE  926 

TYPE  927,  (HEAOER(I),I -65, 96) 

I  J=LKUP( II )-i 

I  C30E=IREADW( 256, HEADER, IJ.ICHAN) 

TYPE  926 

TYPE  927, CHEADERCI) ,1=1,32) 

CALL  CLOSEC(ICHAN) 


1  CONTINUE 
C 

C  ••••*»„„  SELECT  PARAMETERS  FOR  ANALYSIS 

C 

501  CONTINUE 

2  CALL  CLOSEC( I CHAN) 

TYPE  933 

933  ffORNAT( *  TYPE  l  FOR  PROCEED,  0  FOR  D0NE»,3X,t) 

ACCEPT  932 * ISW 
IFdSU.Ea.OGO  TO  999 


954 


30 


9541 


43 


950 


954 


935 


49 

32 


C 

C 

c 


type  934 

CCR3 A  T ( •  SELECT  DATA  FILE,  TYPE  0-9*,3X.S> 

ACCEPT  932  « I SU 
DO  30  1=1.4 
II  =  I3'J*4*I 
NAME  (  I  >  =  INA**E<  I 1  ) 

ICMAN=IGETC( ) 

IL  =  L  OOKUP ( ICHAN. N  AMf  > 

I  COD E=IREADW( 256, HEADER , t .ICHAN  > 

TYPE  911, IDATE.ITIME 
T YPr  9341 

-ORMATC  CHOSE  NO  CF  POINTS, TYPE  9-12*./. 

1  T  l  0  »  *  3 

2  T10. *9 

3  no, *10 

4  TIO.'ll 

5  T10, *12 
ACCEPT  932, IN 
I J  =  2  * • ( IN-1 >  *1 
DO  43  1=1,  IJ 
X( I >  =0. 

Y(I ) =0. 

<=0 

TYPr  950 

-  0  R  V  A  T  ( *  TYPE  1  c  r;o  r  VE  RY  POINT  (DT=.I  SEC), 2  -0® 
1  E  V  r  R  Y  OTHER  POINT  ( D  T  =  .  2  SEC)»»3X,t> 

ACCEPT  932. IDT 
I  S  =  I  D  T  »2 


=  256*./, 

=  512'./. 

=  1024*,/, 

=  2048*,/. 

=  4096*»3X,S> 


I.REC03  =  0 

NRECDS=2**IN/256*I0T 
Ir<NRECDS»9.GT.IL>50  TO  ’97 
TYPE  954 

-OR^AM*  TYPE  1  FOR  HISTOGRAM,  0  FOR  NO',/, 
1*  (ONLY  AVAIL  IF  CHAN  0  SELEC T£ D  )  * , 3X , S  ) 
ACCEPT  93 2  » I HSW 


TYPE  935 

-ORMATd  SELECT  CHANNEL  NO.  TYPE  0-3 • ,/ , T1 0 , • 0  =  3L- 


1 

T  1  0  .  *  1 

= 

ROLL* ,/. 

2 

T  1  0  ,  '  2 

r 

V  STAR* ,/, 

3 

T  l  0  ,  '  3 

= 

H  STAR*  ,/, 

4 

T10,  *4 

r 

V  UNSTAR*,/, 

5 

T  1  0  ,  *  5 

= 

H  UNSTA8*./, 

4 

T10,*6 

r 

COLLINS  ALT*./. 

7 

TIO, *7 

= 

rOLL INS  ERROR  SIG',/, 

9 

T 1  0  ,  *  3 

r 

SPAP-* ,3X,J) 

ACCEPT  93  2 ♦ I C  S  J 
3L0CK  =2 


IF( ICSW.3T.0 ) SO  TO  31 
I  E=0 

DO  43  1=1,22 
'JLCI  >  =  1*200. 

DO  43  J  =1 , 9 
NLIM{ I, J)=0 

IC0DE=IREADU(2304,I NPUT .BLOCK, ICHAN) 
IRECDS=IRECDS*1 

in IRECDS.GT.NRECDS)  GO  TO  1200 
I r < I HSW .NE . 1 ) GO  TO  55 

..........  CALCULATE  HISTOGRAM 
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RADAR •  ,/, 


30  “53  KK  =  1  ,  256  »  IDT 
33  49  1=1.21 

I  F( INQI0(KK).5T.UL( I) )  GO  TO  49 

NLI'MI.1>=NLIM<I»1>M 

30  TO  50 

49  CONTINUE 
NLI*(22,n=Nll*«(22,l)*l 

50  00  52  J  =  1 »  8 
00  51  1  =  1  .21 

I- (  INALG< J.KK) .GT.ULI I > >G0  TO  51 
NLI'*(  I  »  J*  l  )  =  N  L I M{  I.J*1)*1 
30  TO  52 

51  CONTINUE 

I.J*l  >=NLIM<  I»  J*l  JM 

52  CONTINUE 

55  CONTINUE 


C 

C 

55 


33 

C 

C 

c 


31 


34 

1200 


C 

c 

c 


952 

953 
54 

58 

C 


CORRECT  NRL  RAOAR  RANGE 


CONTINUE 

00  33  J  =  1 . 23a  »  15 
<=*<♦  1 

<(<>=IN0I3(J) 

IF(X(<».5T.3000>X(K)=X(K)-4096 
IF(X(K».3T.IOOOJ»IE=IF>1 
JJ=J* IOT 
Y<K>=IN0I3<JJ> 

IC(Y(K>.'3T.3  0  00>Y(K)  =  Y(K>-4096 

I"(Y(<).3T.10000)IE=IE*1 

CONTINUE 


IL0CK=9L0CK*9 
50  TO  32 

I COC E= IRE AOW<  2  30  4  . INPUT. BLOCK, I  CHAN) 
I9EC0S=IREC0S*1 

Ic< IRECOS.GT.NRECDSJOO  TO  1200 
00  34  J= 1,258, IS 
X  =  K  ♦  1 

X(K>=INALO< ICSW, J> 

JJ  =  JMOT 

Y<K) =INALGC ICSW, JJ> 

CONTINUE 
SL0CK=RL0CK»9 
SO  TO  31 
CONT INUE 

IFCICSW.ST.OJGO  TO  1201 
I  Ft  I HSW .NE • 1 ) GO  TO  58 

..........  print  out  histogram 

00  54  1=1,22 

TYPE  952 ,(UL(I>-200.)»UL(I) 

-ORMATC  •  ,  2F8  •  0  ) 

TYPE  953,(NLIM(  I,J),J=1,9> 
r0RMAT(*.*,T25,9I8) 

CONTINUE 

CONTINUE 


& 


¥ 
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c 

c 


PRINT  NRL  RADAR  frrqR  CJJNTERS  "IE* 


TYPE  936. IE 
-0RMAT(»  IE=*  .110  ) 
CONTINUE 


936 
1201 
C 

L* 

C 

937 


938 

37 


946 

C 

C 


33 

947 

C 

n 

c 


36 

951 

C 

r 

C 


939 


38 


C 

C 

949 


PRINT  OUT  RAW  DATA  POINTS 


ty^e  937 

*0R“AT(*  INPUT  NO  OF  POINTS  TO  PRINT*, 5X,S> 

ACC”T  932, NP 

NP=Na/2 

IFINP.EO.OIOO  TO  37 

type  934, < x< 1 > ,y ( 1 > ,1=1 ,npi 

-ORMAT(»  »,10F12.4) 

CONTINUE 
T  =  0. 

TN0=2»»IN 
N0M=2»*{IN-l ) 

TY3-  945, IN, rNO,NCM 

FORMAT!*  IN=*»I4,»  TNO=  *,F6.1»*  NOR:*, 14) 
*»•**.•»»»  CALCULATE  AND  PRINT  «E  AN 


90  35  1=1,  NO'* 

T  =  TMX{  I)  *Y<I 1 1 / TNO 
TYPE  947, T 

-ORRATt*  M  E  A  N  =  ' » c  1 2  •  2  1 
TV=3 .0 


CALCULATE  AND  PRINT  VARIANCE  9ASE0  ON  TI**E  mI^t0pv 


90  55  1=1, NOR 
<m=*u>-T 
Y(I )=YCI)-T 

7 V  =  T  V*X (  I  >»»2/TN0*Y  t  I  )**2/TN0 
TYPE  95 1 , TV  »  T  V  *  •  •  5 

rORM  A  T ( *  0 9  AST D  ON  TIRE  HISTORY,  V 4R I ANCE  =  • , F 1 5 . 3 , 
1  3X,*RI8S=*,F15.3> 


3R INT  OUT  POINTS  WITH  mean  SUBTRACTED 


TYPr  937 
TYPE  939 


rORMAT(  *M2£R0  MEAN>*SX,*> 
ACCEPT  932, N=» 

I F ( NP •FO.OIGO  TO  38 
N3=NP/2 

TY°E  938, ( X( I )  ,Y< I >, 1  =  1 ,NP> 

CONTINUE 

INM=IN-l 


APPLY  R-H  WINDOW.  IF  OESIREO 


TYPE  949 

eORM A  T I •  APPLY  3  TERM  BLACKMAN- HARRIS  WINOOW?  0  FOR  NO. 
I  1  FOR  YES  * . 3X»  $) 

ACCEPT  932, IW 
IF(IW.NE.l»60  TO  45 
■00  4  4  1=1, NOM 
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N=2»!I-i) 

X! I) = 1. 42323-. 49755 *C0S(2*PI»N/TN0> 
1  ♦.07922*COS!2»PI*2*N/TNO))*X! I) 
N=2»!I-l)*l 

V! I)=!.42323-.49755*C0S!2*PI*N/TN0> 
1  ♦.07922*COS!2*PI*2*N/TNO))«Y!I> 


44  CONTINUE 

45  CONTINUE 
C 

C  ..........  CALL  FFT  SUBROUTINE 

C 

C ILL  CWFFT!X,Y ,INM> 

C 

C  ..........  39  I  NT  RAW  'c T  COEFFICIENTS 

c 

Tr°"  o-*o 

740  -CR^ATC  INPUT  NO  OF  RAW  COEFF  TO  °RINT*,5X,i> 


ACCEPT  932, N 

Ir(*J. £3.0)50  TO  37 

30  40  I=l,N 

Ic( I .ST.N0M>00  TO  39 

TYPE  34l,I-l,!I-l>/ITN0».l,ICT),X!I>,Y!I),IX!I)««2.*Y!I>«»2.)...5 


941  -ORMATC  *,I5,4F12.3> 

40  CONTINUE 

39  CONTINUE 

C 

c  ..........  SELECT  NUMBER  FOR  FRE3UENC*  SM00T-IIN3 

C  PRINT  RESULTS 

r 

TYPE  942 

942  -OR«AT(*  INPUT  NO  OF  vCOEFF  PER  FREO  AV5*.5X,*) 
ACCEPT  932  , \- 

TY3E  043 

94  3  rORNAT(»  INPUT  NO  OF  FREO  AVG  TO  °R I  NT • , 5X , t ) 

ACCEPT  932, NPR 
A  T  =  0  » 

00  46  1=2, NON 

AT=ATMX!I)»»2.*Y!l ) *  *2 • > /2. 

45  CONTINUE 


00  41  1  =  1, 8000, NF 
AA  =  0  . 

00  42  1  1  =  1  , NF 
IF(  IMI.GT.NCMJGO  to  9q  7 
AAiAAMX!  I*Il  )**2.*Y!I*II>**2.)/2 
42  CONTINUE 


995 

47 

944 

41 
99  7 
998 


TYPE  944,1  ,  1 ♦  N-- lit  I* ( NF- 1 ) / 2 . ) / ( TNO • . 1  *  I OT ) , AA 
IF< ICS W.LT .2)50  TO  47 
I  Ft  I CSW • 5T .5)50  TO  47 
TRL=<TN3*.W5T)/(  I*(NF-l)/2.) 

AMC=!TRL/!2.»PI) )  »•  4 

AO=A A»AMF 

TYPE  995, AO 

FORMAT! •  ♦  • ,T  43, F 15*4) 

CONTINUE 

FORMAT!*  *,2I5,F12.3,3X,F15.5) 

IF1I/NF.3T.NPR)G0  TO  997 
CONTINUE 
TYPE  998 

FORMAT!*  NOT  ENOUGH  POINTS*) 

TYPE  996, AT ,AT  »*.5 
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•3* 


& 


& 


<r 


& 


o  o  o 


936  -ORMATCOTOTAL  V  AR  I  ANCE=  •  ,  F 1 5.3 , 3X  »  '  ■<  M  5=  •  ,  FI  5 . 3  ,/// ) 

•»*••*••••  LOOP  TO  BEGINNING  OF  PROGRAM 

30  TO  2 
999  CALL  CLOSEC(ICHAN) 

CALL  I  FR E E C( I  CHAN ) 

STOP  'ALL  DONE* 

EMO 
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noonoooo 


=»sr  FOURIER  TRANSFORM  ALGORITHM 

(REF;  RICE,  PE  TER  K.t  "SPEED  UP  THE  FFT". 

ELECTRONICS  DESIGN,  APRIL  2S.1970) 
ADAPTED  9Y:  DAVID  A.  WALDEN 


4 


7 

9 


II 


SJHOUTINE  DUrFT(X, Y,IS) 
INTEGER  O.Gl 
DIMENSION  X(l) »T(1 ) 

3=1  o 

N  =  2  *  *  3 

=  =<• ATANI 1 . >/N 
DO  4  1=1, N 
X(I>=X(I)/N 
rm=Yd)/N 
CONTINUE 
00  9  LI  =  1  , G 
L  =  L  I  -  1 

31  =  2  * • ( 3-L - 1  ) 


L  2  =  2 • *L 
00  9  1  =  1, L2 
U=v/Jl 
<1=INT(R1  ) 

<2  =  0 

00  i  <=1,3 
4  2  =  <  1  /  2 

<3=<1-2«INT (22) 

<  1  =  I  N  T  (  p  2  ) 
t-(<i.E0.3)30  T0  ' 

<  2  =  <  2*2* • (3-<  ) 

C3NT  INUE 

Y l =C  OS (P*  <2 ) 

Y2=-3IN(P»K2) 

00  7  J  =  l  ,  51 

Y3=Y(M*31*1)»y1“Y(M*G1*1)*Y2 
Y4=X (M*5l*l) •Y2*Y(M4Gl*l)*Yl 

X(M*G1*1)=X(m*1)-Y3 
Y(M*31*l)=Y<M»l)-YA 
X(M*1 >=X(M»1>  *Y3 
Y(M*1)=Y(M*1)*Y4 
M  =  M*  1 
CONTINUE 
M=M*51 
CONTINUE 
00  10  L3  =  l ,N 
I =  L  3  —  1 
<1  =  1 
<2  =  0 

00  11  <=1,G 
R4=< 1/2 

<3=K1-2*INT(R4> 

<1 = I NTI R4) 

IF <K3. E3.0>50  TO  11 
<2=<2*2**«G-K) 

CONTINUE 

IF(<2.SE.I )GO  TO  ip 


& 


.  , ' 
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*K3=X< 1*1 ) 

X(  I*U=X(«2*1> 

«<K2*1)=RK3 

*K3=y<I*l> 

r<  I*l)=Y<K2*l) 
y(K2*n=RK3 
iNuc 
L  3-N  /  2*  l 
TO  13  1  =  1, L5 

<-r-  i 

J  =  I-2 

«  =  X(  I  >*X(N-J) 

^  =  (Y(  [I*Y(N-J))  •r0SCK»P/2)-(X  (I  )  *X ( Y-J )  >  *SI\(K*?/2) 
-  =  Y<  I  )-Y(N-J) 

)=(Y(t»*Y(N-J)»«SIN(K*P/?)«'<X(I)-X(N-J>>*COS(K»P/'») 
X( I)=(ft*5)/2 
Y( I >  =(C-T»/2 
XCM-J)=( A-R)/2 
Y ( M- J ) = { -C -T ) / ? 

COMTI.Yjr 
* -tlmn 
-:no 


F/G  1/8 


AQ-A0G3  6H7 
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